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SYNOP5E5 


of the 

Ph- D- Dissertation 
on 

A SPATIAL INTERACTION MODEL FOR STATISTICAL ANALYSIS 
OP SOME INTERSPECIFIC COMPETITION EXPERIMENTS 

■by 

M. Narayana Reddy 
Department of Mathematics 
Indian Institute of Tecnnology 
Kanpa r— 20801 6 (India; 

Interspecific competition experiments involve growing 
two or more species togetiner* 'When these species are grown 
in separate rows within a plot then it is called intercropping (ic). 
Among interspecific competition studies intercropping is one 
of the important topics (Mead/ 1979) . Currently there is a 
growing interest in IC studies because of its substantial 
advantages over mono croppincf (mc) / i.e» growing them In separate 
plots* The importance and research needs of IC were reviewed by 
Willey (1979) # while Mead and Riley (1981) reviewed the existing 
methods of statistical analysis* The conventional methods of 
statistical analysis associated with IC are not satisfactory due 
to the complex nature of interactions involved* Most of the 
methods suggested so far emphasize on testing the difference 
among the treatment means- However# for studying the effects of 
some important factors# sucb as spatial row arrangements and plant 
densities # fitting a relationship py incorporating these factors 
is important. This helps ±n arriving at the optimum row arrangement 
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and plant densities* In the present study efforts have been 
made for developing suitable relationships for two component 
crop experiments involving various row arrangements and plant 
densities of both the crops* In the first chapter motivation 
for the present study is discussea ana me relevant literature 
is reviewed. 

Different row arrangements in IC introduce different degrees 
of spatial interaction among the plants or the crops* It is 
of interest to find the row arrangement which maximizes a 
meaningful linear combination, of two outputs* For this relationships 
are formulated for individual row yields based on the principles 
of plant competition* When rj^ rows of one crop is followed by 
r 2 rows of the other crop t for fixed intra and Interrow 
distances • the following relationships are assumed - 

ni “ ”1 + hk yk,i-i + hm + “li ' i = a - 

k = 2 when i = 1 / m = 2 when i = r^ » otherwise k = m = 1 

and = ^^k-r^+rj ' 

Yjj = ^2 + S 2 J, *• + Ujj . j = : 

k =: 1 when j = r^+l# m = 1 when j = ' otherwise k = m = 2 

ym.r^+r^+l = ' 

In the above equations » k = 1/2/ is the expected 
row yield of crop k in the absence of competition '• y^^^ and 
y 2 j are the random variables corresponding to the yield of 
the first and second crops / respectively / 


is the 



xi 


competition coefficient between any two nearest neighbour (NN) 
rows of crop k * ' k ^ m = 1#2 represents the competition 

of crop m on crop k t and U 2 j are errors and assumed to 

be independently distributed as N(0#a^} and nCO/CT^^^ respectively. 

The above equations belong to the simultaneous spatial 
process introduced by Whittle (1954)- A bivariate model is 
obtained for the plot yields from the above set of equations- 
The parameters in this model are functions of row arrangement 
and proportion of the crops- Models for MC situations in terms 
of their plot yields are also developed- Some alternative 
formulations » such as conditional autoregressive model 
(Besag# 1974) # are also discussed- In the present situation 
simultaneous model appears to be more appropriate and hence 
it has been considered for further studies- The development 
of models for row arrangements is discussed in the second 
chapter- 

The problm of estimation of parameoers has been discussed 
in the third chapter- For this purpose the method of maximum 
likelihood (ML) is adopted in estimating the parameters* 
Estimators for replicated field experiments are obtained using 
the information from IC and MC plots- As the model is nonlinear 
iterative equations based on Newton-Raphson method for obtaining 
the ML estimates have been developed. Expressions for obtaining 
the large sample variance-covariance matrix are derived. The 
model has been extended to take care of block effects- To 
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examine the performance of the ML estimators * simulation study 
has been carried out by assuming some practical experimental 
situations* The effects of changing the error variances and 
the size of the experiment have also been studied* 

In addition to row arrangements the magnitude of competition 
depends upon the level of plant densities* In the fourth chapter 
the model is extended to incorporate the parameters that quantify 
the effects due to plant densities* For this* the relationship 
between competition coefficient* say P * and the NN interplant 
distance x * is assumed as 

B(x) = , 

where and c are constants*A mernoa for, obtaining the ML 
estimators of the parameters is developed* Expressions for 
large sample variance - covariance matrix are also derived * 

The results of a simulation study are also discussed* 

In the fifth chapter some methods for obtaining the 
optimum row arrangement and plant densities are discussed * 

Certain linear combination of output such as land equivalent 
ratio (LER) * in terms of expected yields » are considered 
for optimization purposes * For obtaining the optimum row 
arrangement * comparison of expected output for various row 
arrangements * appears to be a convenient method * However * 
for obtaining optimum densities nonlinear optimization metixjds 
have to be used* These methods are illustrated with some examples* 
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In the sixth chapter the model has been fitted to an 
experimental data. Certain aspects of the adequacy of the 
fitted model to these data have been examined. The overall 
fit of the model appears to be satisfactory- The advantages 
of the fitted model / in gaining some insight in the nature 
of competition » and optimum row arrangement and plant 
densities / have been discussed- 
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CHAPTER 1 


INTRODUCTION REVIEW AND SCOPE 

. 1 INTRODUCTION 

In plant ecology competition occurs among plants when they are 
onsuming limited resources* In such situations the presence of any 
lant changes the environment of its neighbouring plants which may 
ffect their growth rate and form- Different forms of competition 
xperiment and some of the related biometrical problems have been 
eviewed by Mead (1979) - The purpose of competition studies is to 
nvestigate the competition relationships among plants of the same 
pecies or that of different species i-e- how the plants modify their 
ields due to joint utilization of the resources when they are g 3 X>wn 
ogether- These relationships add greatly in understanding crop 
rowth and production* 

An interspecific competition study is the one in which plants 
f two or more species are grown together on the same plot* In such 
tudies usually yield of each species is recorded separately* The 
Dst complex and important among the interspecific competition studies 
s intercropping* Growing two or more species together in separate 
3WS in which various arrangements of rows are possible# is called 
itercroppinq (IC)- The importance of intercropping# its agronomic 
incepts and research needed along with some concepts in the assessment 
£ yield advantages* have been reviewed and discussed by Willey (1979). 
le study herein is aimed at developing suitable statistical methodology 
3 r analysing the effects due to some of the important factors in 
itercropping involving two species* 
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More generally, growing crops by irregular broadcasting or 
by mixing witnin the rows is called mixed erpppinq- Mixed 
cropping is an age-old cropping system in India and also in 
the semi-arid tropics throughout the world { Jodha# 1979) . 

The subsistence farmers in these areas practise mixed cropping 
because of its significant advantages* Some of them are the 
following • 

1* It reduces the risk of total failure of crops per unit 
area in unfavourable weather conditions* 

2* It utilizes the limited environmental resources like land/ 
water and nutrients etc*# more effeciently and may result 
in larger output than growin^a tnam separately on the same 
piece of land. 

Intercroping is a modification of mixed cropping where each 
component crop is gro'/ai in separate rows. 

Component crop 155 either of the individual crops making up 
the IC ) intercrop yield is the yield of the component crop in 
IC expressed over the total intercropped area (i*e* area occupied 
by both the crops ) ) Monocropplnq (MC) or sole cropping refers 
to a component crop being grown alone at optimum plant population 
and spacing / unless indicated otherwise* Some examples of IC 
and MC are sho''/m in Figure 1*1* 

The interrow (between rows) and intrarow (between plants within 
a row) distances for any component crop in IC need not be the same 
as that of MC. When an IC treatment is obtained by replacing certain 
proportion of a monocrop with an equivalent proportion of the other 
crop / it is termed as a replacement series treatment* Replacement 



0^ 

H 

H* 

p 




h 

d 

d 

d 




rp 


0 

a 




H 

M 


P- 



o 

d 

g 


vQ 


(t 


pb 

o 

t2S 

• 



ft 

p 

d 

0 


X 


cr 

CD 

0 

H 

H 

CD 

CD 


0 

3 

« 

p 


o 

P 


P 

H 

d 

V. 


<3 

H» 


Pu 

0 

0 

H 

tJ 

H 

C/) 

(A 

§ 

P 

(A 

K 

0 

O 

H* 

3 

H 



3 


ft 

0 

ud 

0 

ft* 

P 

P 

H* 

d 

P 

Hi 

H 

d 

0 



3 

P 

P 

CD 

t3 

(t 

d 

Pu 

X 




c 

P 

H* 

ft 

0 

(D 

3 

d 

ft 


d* 


P 

u 


H* 


p 


d* 

P 


0 

H 

it 

3 

d 

CD 

d 

P 



CD 



P 

CD 


CD 

d 

O 

P 

M 

p 


ft 

hh 

O 


0 

d 

b 

C/3 


H* 

0 

Hi 

ft 

O 

P 

ft 

P 

0 


CD 

t 

d 

W 

3 

P’ 


(A 

h 

P 



d 

0 

w 

p 

P 

o 

P 

ft 

hh 

P 


cr 

H 

ft 

P 



3 


P' 

H 


d 

p 

P 

iQ 

ft 

ri 

a 


d 

P 

p* 


(3 

V 

H- 

ft 

d 

0 

0 


O 

w 

0 

d 


p 

H 

d 


ft 


V 

p‘ 

P 

w 


a 

P 

p- 

h. 

d 


0 


d 

H, 

rt 

CD 

d 

P 


iQ 

P 


"d 

P 

p 


h 

P. 

P 

P 

% 

H 


p 

CD 

ft 

iQ 


P* 

H 

d 

d 

K 


3 

O 

ft 

03 

p 

P* 

d 

P 



H- 


< 

p 

d 


CD 

ft 


P 

H 

ft 

P 


H- 

Q 

d 

P 

(D 

d 

P 

CD 

0 

0 

P 


Pj 

o 

W 


W 

P- 

K 



d 


d 

3 

P 

0 

P 

0 

*d 

H 

0 

0 

CD 

Hi 

h 

O 

M 

d 

• 


h 

0 


d 

0 


ft 

P 

P 

0. 

o 


tf 

d 

0 


p 

H 


P 

ip 

3 


0 



tf 

P 

3 



0 

§ 

H- 

H 

P- 



0 

p 

d 

t 

X 

p- 


3 

d 

P 

d 

<+ 

d 



ft 

ft 

2 

d 

iQ 


0 

w 


p 

H 



d 


0 

rt 

P 



p 


d 


CD 

S 


d 





O 


ft 





« 



X- -X — X--X--X 
X — X — X — X — X - 

s ^ ^ — x- -X- “X -- -X - 

O Pj 

^ X- “X - “X - -X- -X - 

X- -X - -X — ^x--x- 
X- -X--X~-X — X- 

o- -- o — o — o — o- 

O 0~ ~0 — O O' 

X ^ 0 -- 0 — 0 -- 0 -- 0 - 

Q (D 

M o--0-~~0 — 0--0- 

0--0--0--0- -o - 

C >__0 — O — O --0 - 






X- 

-X- 

-X- 

-X- 

-X 

--X 

tf 

t1 

H- 

0/ 











(P 

P 

P 












ft 

n> 




o- 

--o - 

-o 

~o - 

“O 


cf 

H 

(D 












cu 

h 

H* 










(D 

0 

0 

3 




x~ 

-X- 

-X- 

-x- 

-X 

--x 


(t» 












tn 

3 

0 

0 










Q) 

(D 

t:i 

p 




o- 


-o * 

-o- 

-o 


B 


0) 

rt 










(t) 

ft 


03 












0 

P 




X- 

-X- 

-X-- 

-X- 

-X 

X 

ti 

03 

< 

ct 










0 

(D 

(D 












h 


H- 




o- 

-o- 

-o~ 

-o- 

*’ o- 


& 

H* 


D 











(D 

3 











h 

0) 

0 

P 












P 

P 










0 

ft 

0 

a 




x- 

-X" 

“X- 

™X“ 

-X 

- -X 


h 

0 

(t) 











0) 


»i 










n 

OJ 

0 

03 




X" 

-X- 

“X- 

-X- 

~x 

--X 

d 

ct 


ft 










t 

B 

03 

a* 










03 

0) 

• 

p 


H 


o- 

-o- 

-o- 

--0- 

-o 

—o 




Oj 


o 

cr 







0 

ft 



— o 









»-h 



p 




X- 

•-X ' 

--x- 

-X- 

-X 

--X 



H- 

iD 










ft 


iQ 


— o 









If 

(D 

P 

p3 




X- 

-X- 

-X- 

-X- 

-X 


(D 



p 











(D 

(D 

Of 

“ -O 









0 







o- 

-o- 

_o- 

-o- 

-o 

~-o 

rf 

0 


CD 










D* 


# 

< 

— o 









n) 

0 

H 

03 











'd 

lO 

H 

r* 

— o 



X- 


~~Jxr“ 

-x- 

~x 

~ -X 

o 


H- 

M 

0/ 










M 

0 

< 

rt 










o 

>: 

(D 

H» 

--0 



X- 

-X- 

-X- 

-x- 

”X 

-"X 


01 

C/D 

P 










« 



<Q 











H- 

(1/ 





a-- 

~0“ 

-o- 

-o~- 

-o 

- -o 


d 

P 




o 










H- 



V--' 








DJ 

CD 

CD 




a- 

~0" 

-o- 

“O— 

-o 

-o 




H 











:s 

s 

a 

Cb 




x~- 

-x™ 

-X-* 

-X-- 

-X” 



\ j 


03 












H 

Qi 











h 

(t) 

<j 




X- 

-X- 

-X-- 

-X-- 

-X' 

--X 


(D 


03 












O 

P 











H 

Hx 

ft 











(D 


03 












U) 

iQ 











H* 


fD 











0^ 

03 

0) 











0 













CD 

P 

0 











a 

















w 





4 


:rops* Two chief objectives of these experiments (Mead and Willey^ 1980 
*earce and Gilliver/ 1978 / Vandermeer# 1986^ are the following • 

To determine v.hether a given IC combination is indeed better than 
MC- 

To understand the underlying mechanism in IC for enabling further 
improvement- 

To meet these objectives a research worker is required to 
level op appropriate models for analysing data from experiments 
.nvolving spatial row arrangements and plant densities- The present 
;tudy concentrates on these aspects or the IC experiments* 

.2 REVIEW OF THE LITERATURE 

.2.1 Existing methods of statistical analysis 

A good account on the current raetnods of statistical analysis 
.f IC is discussed by Mead and Riley (198i; and these are further 
llustrated by Dear and Mead (1983<’1984) - The statistical methods 
rhich are in current use are mostly for testing significant differences 
mong the treatments by analysis of variance (ANOVA) or covariance 
ANCOVA) techniques- Univariate methods have been mostly used for 
nalysing the individual crop yields and some derived variables such 
s monetary value and nutritive value- Other important variables 
fhich have been often analysed are some indices of the combined 
’iel ds which characterise the competition. The most popular among 
he competition indices of the combined yields is 'Land Equivalent 
:atio' (LER)/ as discussed by Mead and Willey (1980) • LER is defined 
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as (Willey and Oslru-1972J 

LER = + L 2 = ^ ^ 2 /^ 20 ' (1.2.1) 

where Y^,/ (ic = 1/2) is the yield of the component crop h in IC/ 

is its corresponding MC yield and ip the partial LER of crop k, 
i.e. proportion of the kth component crop yield to its raonocrop yield- 
LER is widely used by agronomists because of its interpretation as 
the relative area required to produce the same yields that are 
achieved in IC through MC • In the LER analysis one of the important 
hypothes.es to be tested is * LER = i against : ler^ 1 , and 
this is discussed by Skovgaard (1986). When is accepted the 
yield advantages it IC are not significantly different from MC- 

In IC the growth rate of any component crop is not independent 
of the other. Hence analysis of ic is essentially a bivariate 
problem/ as analysis by univariate methods may not detect the real 
treatment differences- Pearce and GiiHver (1978) suggested a 
standard multivariate technique of analvsis of dispersion for testing 
the treatment effects- They have demonstrated the use of graphical 
procedures for comparing treatment mt^oiis for the two species 
simultaneously- These graphical procedures have been extended 
(Pearce and Gilliver- 1979 / Gilliver and Pearce, 1983) (1 ) to compare 
the advantages of IC over MC, (2) to identify the IC system that 
naximizes the total produce when there are certain constraints on 
the component crop yields/ and ( 3 ) to study the two and three factor 
Interactions in factorial experiments- The analysis suggested by 
Pearce aid Gilliver is based on the assumption that the correlation 
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coefficient remains the same for all the trecitrnents which may not 
be true in many situations* To overcome this Singh and Gil liver (1988) 
suggested an extension to the analysis suggested by Pearce and Gilliver 
by considering the corrrelated error structure among the treatments* 

When the IC treatments Include certain factors such as various spatial 
row arrangements and plant densities of the component crops/ it is 
important to develop appropriate response relations among the component 
crops yield rather than looking at the significance of certain 
treatmental differences* Such relationships will be of great help 
in understanding the response of yield in IC to varying spatial row 
arrangements and plant densities# This relationship may be useful 
in finding the optimum IC t 3 reatment* 


Vandermeer (1986) suggested a relationship in terms of plant 
yields of component crops to predict the yields for different spatial 
arrangements of the plants in IC. The suggested equation is 


w. . ssA-j •- S 8^. .m, . "*2 8 j 

11 Oli jtl ^lij Ij / 1 ^2ij 2j 


( 1 . 2 * 2 ) 


where is the yield of the ith plant of the first crop v 

is the expected yield attainable in the absence of competition * 


^lij represents the competition effect of the jth plant of the first 
crop on the ith plant of the same crop / 

competition effect of the jth plant of the second crop on the ith 


plant of the first crop/ and 0 ) 2 j yield of the jth plant of 

the second crop* Some adhoc methods have been used to estimate the 
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parameters in the above equation through the individual plant 
data from small systematic designs* In IC the total yield of 
a plot (or the average plant yield) is important# hence the 
statistical analysis based on the model that explains the 
variation among plot yields of the component crops# is more 
realistic. 

Regarding the yield density relationship "Wright (1981) suggested 
the following models for analysing the experimental data of crop 
mixtures# 

^1 “ ^ll'^ + ^12^^ ^ ^1 (1.2.3a) 

^2 * ^2 ^ “2 ^22^2 ^21*\^ ^ ^2 (1.2. 3b) 

vjhere is the yield per unit area of crop k# and are the 

plant densities of crop 1 and crop 2# respectively # (lc#ni = 1#2) 

rn 

are the competition parameters# and are the errors* Here 
Y-^ and Y 2 are independent 'vdiich may not be true in IC. 

Another important aspect of IC analysis is to compare the 
stability of IC treatments on crop yields and monetary value over 
years. This aspect is not considered in the present study. Some 
of the papers dealing with this analysis are by Pearce and Edmondson 
(1982#1984)# Mead etal* (1966) and Singh etal. (1988) . The designs 
that are currently in use are randomized block design (RBD)# split 
plot and systematic designs (Mead and Stem# 1980) . 

As the paresent study is concerned with the statistical modelling 
and analysis of IC experimental data# the methods for analysing 



raditional competition experiments like crop mixtures are not 
evlewed here. Some of the papers in this area are by Williams (1962) . 
cGilchrist (1965)/ KcGilchrist and Trenbath (1971)/ Federer (1979) 
nd Federer etal. (1976). 

.2*2 Spatial interaction models 

.2.2*1 Introduction 

When a set of spatially located variables / i =l/2/.*./n/ 
iisplay interdependence over space then the data are spatially 
I utocor related. The autocorrelation may be due to reaction or 
Interaction among the variables* Adjacent: trees may compete for reso- 
urces like water /nutrients and suniiynt etc-/ displaying between tree 
Interactive effects/ but they may also react to the general 
availability of nutrients within the reach of their root system* 

When reaction is dominant/ a regression model is appropriate/ 
whereas interactive effects suggest the need for a model with a 
spatially dependent covariance structure. The model should include 
the parameters of reactive and interactive effects as both are 
important factors* A.s the present study is concerned with developing 
the statistical models for IC / based on spatial interaction among 
the plants/ the spatial interaction models are briefly reviewed 
in the following section. 

1. 2*2*2 Simultaneous model 

Whittle (1954) proposed a simultaneous spatial autoregressive 
scheme of the following form 

X. = P. -f E s. .(X.-p..) + £/ i/ j = 1/ 2/ . . / n ) (1.2.4) 

j;^i J J 

where E(t .) = o and cov(t ./£ .) = oj ir r = j/ zero otherwise. 

1 i J i 



the matrix notation this can be expressed as 


X=4+S( X-M 

lere X = (x,#***# X 

— 1 n 

a = iz ^ f * * * t E * 

— 1 n 


4- e # 


M = ( / • • • # ^ 

^ = ‘^ij'nxn =il 


5 then have 


0 . 


(1. 2.5) 


E(£) =0# V(£) = A = diag (cf? a^) . 

— — ^ 1 n 

E(X) = M # Vg(x) = (I-S)'"^ A(I-S')“^. (1.2.6) 

ne necessary and sufficient condition zor the existence of 
his process is (l-s) should be nonsingular. In the present study 
t is assumed that X is distributed as multivariate normal with 
ean 4 and variance V^* 

• 2.2.3 Conditional __ model 


This approach requires to specify the conditional distribution 
>f each random variable X. given X. = x. • j i. In the case 
>f nearest neighbour situation this approach was first suggested 
3y Bartlett (1955) and the corresponding mathematical and statistical 
theory was developed by Besag (1974.1975) • Suppose that the 
zonditlonal distribution of X^ is normal with conditional mean 


E(XJX. = X. . j i) = 4. + S c.. (x. - M-,.) 
1 J J 1 J J 


(1.2.7) 


and conditional variance 


Var (X^ I Xj 


X . . 

3 


j i) = aj 
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where = 0. then the dispersion matrix of the scheme is 

given by 

= (I-C) ^ A , where C = (c.J . 

ij nxn IX- 2. a; 

The scheme exists only where V^is symmetric and positive definite- 
Besag (1974) named the above model as auto-normal model- In fact 
he introduced a general class of auto-models based on the conditional 
probability formulation- Here X has a multivariate normal 
dist:ribut:ion with mean 4 and variance v * 

— Q 

2* 2* 2* 4 Difference between simultaneous^ and conditional models 

The general distinction betv.een these two approaches Is due to 

Brook (1964). The difference lies in the specification of the 

dispersion matrix- In the case of Gaussian process Besag (1974) 

discussed a simple method of deriving an equivalent conditional 

scheme from simultaneous scheme. Bartlett (1975) discussed the 

same through spectral density functions- In the case of Gaussian 

2 

process v/hen A = o l . the simultaneous and conditional schemes 
will be identical only if (Riply/1981j 

(I - C) = (I - s') (I - s) - 

In the case of simultaneous model there is no unique representation 

of dispersion matrix in terms of S/ vfiereas in the conditional 

case the dispersion matrix determines C uniquely - But the 

major disadvantage vjith conditional anproach is about their 

—1 

existence because of the condition that A (l-c) should be 
symmetric and positive definite- These rwo models are named as 
spatial autoregressive models analoguous to those in time series- 



11 


1*2*2*5 Moving average models 

As in -the case of time series (3ox and Jenkins# 1970) moving 
average (MA) models can be specified by taking = O without 

any loss of generality# as follows (Cliff and Ord# 1981)# 

X. = + 2 a a , i,j = l#2#-..#n : (1*2.9) 

^ ^ j/i ^ 

where E(£^) =0# Var(£^) = and Cov(£.^#£^. ) = 0 where i ^ j- 

In the matrix form these can be e:cpre^=ed as 

X = (I+A)£ # where A = # 

Var(x) = (I-p^) A (I+a" ; . (1*2*10 

The process is representable by simultaneous spatial process only 
if the largest eigen value of A is less than unity in absolute value- 
Haining (1973a) has discussed in detail about these schemes- MA 
models have been mentioned often in literature but rarely used 
in practice (Haining# 1973b) in spatial studies- 

1-2*3 some applications of spatial int eraction models and methods 
of estimating parameters 

1-2* 3*1 Applications in field experiments 

Spatially autocorrelated data occur in a wide variety or 
scientific disciplines# especially in ecology and geography- A 
good description on analysing these type of data is given in the 
books by Bartlett (1975)# Riply (1981)# Cliff and Ord (1931) and 
Upton and Fingleton (1985)* Recently there has been a considerable 
interest in the applications of spatial autoregressive models in 
agricultural field trials- Whittle (1954) and Besag (1974#1977a) 
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fitted these models to analyse uniformity trial data. This type 
of analysis helps agronomists in planning their future experiments- 
Making use of autoregressive models Bartlett (1978) reexamined 
theoretically the method suggested by Papadakis (l937j/ for 
adjusting the treatment means based on the neighbouring plot 
residuals* 

Mead (1966.1967#1968/1971} adopted simultaneous approach 
to study the plant competitions in MC through small competition 
experiments* He fitted the following model for the data from 
systematic design with hexagonal arrangement of plants (Mead/1967) - 

w. = 5 + 2 (w.— 5)-j-e. . (l«2«ll) 

^ jCi) ^ 

where j( i) is the set of neighbours of the ith plant# is the ■ 
yield of the ith plant# 5 is the population mean# 6 Is the 

# 

competition coefficient among the neighbouring plants and s 

are errors v^ich are independently distributed with zero mean and 
variance 

Kempton (1982) fitted simultaneous model for incorporating 
the competition effects due to neighbouring plots in the variety 
trials* The model is 

Y. = T + P 2 Y./p + £. # (1.2*12) 

ir r j xr 

where j(i) is the set of p neighbours of plot i# Y^^ is the yield 
of the rth variety in the ith plot# P is the common competition 
coefficient among the neighbouring plots # and errors 
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independently and normally distributed with zero mean and variance 
2 

d . The above model with block parameter*? to account for large 
scale fertility differences has been fitted by 3esag and Kempton (1986) 

1*2*3*2 Estimation procedures 

In the case of simultaneous models ordinary least squares 
estimators are inconsistent as pointed by Whittle (1954)* Oird (1975) 
suggested a modified least squares method which yields consistent 
but inefficient estimators* ^nhereas in conditional models least 
squares estimators are consistent* Besag (1972*1974*1975) suggested 
coding and pseudo-likelihood (Besag* 1975^ 1977b) methods for 
conditional models* In the normal cases pseudo-likelihood estimators 
correspond to the least squares estimators* In the case of 
maximum likelihood method* the log-likelihood for simultaneous 
and conditional models represented bv the equations (1*2*5) and 
(1*2*7)* when (i = l*2*...*n) is given by 

21og L = -n. log (2cr^) + loglBl - o^(X-A ) ^ B(X-i(£)* (1*2*13) 

where B = (I-S^j (I-S) for simultaneous Ihodel and B = (l-C) for 
conditional model* and 131 is the determinant of 3* 

If h- is of the form D0 * where D is the design matrix* 
the maximum likelihood estimators of 9 * u and the parameters 
in B * are given by 

^ = (d'BD)”^ d'b X (1*2*14) 

n'‘^(X-DS) B(X-Dji) * (1*2*15) 

The estimates of parameters in B may be obtained by minimizing 
-n"^ log IBl + log { n"^ (X-DS) 'b(X-D9) } * 


(1*2*16) 
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Estimation of 0 and B proceeds iteratively by successive use of 
equation (1*2*14) and minimization of {l.2.16i. Computationally 
the problon arises in the evaluation ox iBi at each stage when 
minimizing (l.2*16J iteratively. When B involves only one parameter 
e.g. 5 = / where W is the matrix with bnovm coefficients / Ord (197e 

evaluated IBl as 

n 

lEi = 7T (1-3X. ) ^ 

i = l ^ 

where s are the eigen values of VJ. Kempton (1982) also used 
this method in evaluating IBl. Mead (1967) in small sample 
competition studies evaluated IBl by Lanczo^s (1957) method of 
interpolation which is based on Chebyshev polynomials. In the 
case of simple simultaneous models involving one or two parameters 
Whittle (1954) used general Fourier inversion method for evaluating 
IBl when large number of observations are available on lattice 
systems. 

1.3 SCOPE O? Tim PRESENT STUDY 

When the resources are uniformly dic^tributed then the yield 
variations of the component crops in IC are mainly due to the 
changes in the degree of plant interactions among the plants 
of the crops and random errors- Hence the statistical analysis 
based on spatial interaction models which incorporate the spatial 
dependence into their covariance structure# reviewed in section 1*2.2# 
appears to be more suitable for analysing IC ejqperimental data. 
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In the second chapter model incorporating effects due to 
various row arrangements of the component crops# based on the 
principles of plant competition is developed. Model is formulated 
in tepas of individual row yields. Model for plot yields is derived 
from the model formulated in terms of row yields since plot yields 
are usually observed in IC experiments* Model corresponding 
to MC is also derived* In section 2*3 some alternative models are 
formulated. Similtaneous model appears to be more suitable for 
this situation* Hence this model is considered in details in the 
subsequent chapters* 

Estimation of the parameters through a replicated field 
experiment is discussed in the third chapter for the simultaneous 
model. The method of Maximum likelihood (l-lL) is used for estimating 
the parameters* The ML estimators are derived in section 3*2 using 
information from IC and MC plots* As the model is nonlinear iterating 
equations based on Newton-Raphson method have been developed in 
this section. Expressions for asymptotic variance-covariances of 
the estimators are derived in section 3.3» The model is extended 
to incorporate the block effects in section 3-4 and cor ire spending 
estimators are derived* The performc-nce of ML estimators are 
examined through simulation study and the results of this study are 
discussed in section 3*5. The effects of change in error variances 
and size of the experiment are also studied in this section* 

In the fourth chapter the model is extended to incorporate 
the parameters that auantify the effects due to plant densities* 
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This extension is discussed in section 4* 2* Model corresponding 
to MC is also developed* ML estimators for this extended model are 
derived in section 4* 3* In section^* 4 hlock effects are incojrpo rated 
into the model. The results of the simulation study carried out for 
examining the performance of the i-’iL esriiuacors are discussed in 
section 4.5. A particular case of the model for variation only 
due to plant densities is discussed in section 4* 6. A numericcd 
method for obtaining the ML estimates by iteration and the ej^ressions 
for the asymptotic variance-covariances of the estimators are given 
in the appendix of this chapter* 

In the fifth chapter some methods for obtaining the optimum row 
arrangement and plant density levels# using the models developed 
in earlier chapters# are discussed. These aspects are discussed 
in section 5-2. ' Some examples are given illustrating the method. 

In the sixth chapter the model developed in the fourth chapter 
has been fitted to some experimental data. Some aspects of the 
goodness of fit of the model are also examined* The interpretation 
of the fitted competition coefficients and their use in giving 
some indication regarding the optimum row arrangement and density 
levels are also discussed. 
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1.4 NOTATION AND ABBREVIATIONS 

ML maximum likelihood 

IC intercropping 

MC monocropping 

NN nearest neighbour 

LER land equivalent ratio 

LERg land equivalent ratio in expected yields 

MER monetary equivalent ratio 

MER monetary ecuivalenr ratio in expected yields 

e 

L, partial lend equivalent ratio of kth crop 

iC 

L-. partial land equivalent rario of kth crop in expected 

yields 

total monetary value 


TMV 



CHAPTER II 


SPATIAL INTERACTION MODELS FOR CROPS VJITH DIFFERENT 
I'DW AKR<-\!4GEME:-)Ti 

2-1 INTRODUCTION 

The eiavantage of IC over MC is due to tne cidj us -ability 
of the component crops in a given environment such that the 
total competition pressure per unit area is reduced- This 
suggests building up a model for IC systems in terms of 
parameters which quantify the competition effects in different 
situations- The reduction in total competition pressure in 
IC is due to less intensity of interspecific competition 
than intraspecific competition (Trenbath/19'7 4J Willey# 1979) - 
In other vx)rds the intensity of plant interactions between 
species is less than that of within species- The arrangement 
of the plants on the plot is a basic factor that modifies 
total competition pressure- Hence spatial arrangement of the 
crop rows is one of the major factors that is being studied 
in IC research programme currently- Some appropriate spatial 
interaction models in terms of the competition parameters for 
the row arrangement of the component crops in IC have been 
developed and discussed in this chapter- 

2-2 DEVELOPMEi'IT OF SIMULTANEOUS MODEL 

2-2-1 First-order simultaneous autoregressive model 

Let us consider a plot with rows of crop 1 and rows 
of crop 2- The arrangement of these N = (N 2 ^+N 2 ) rows are 
such that r^ rows of crop 1 are rollowed by r 2 rows of crop 2 
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for resources is likely to be severe between nearest 
neighbouring (nn) rows* In developing the nK>del the following 
assumptions are made. 

1. The competition is limited only to the NN rows and its 
effect is linear on the NN row yields- 

2. The first and the last rows are neighbours* 

To implement the last condition we assume that appropriate 
guard rows are taken in field e:<periment. 

Following these assumptions the individual row yields 
in an of the component crops can be expressed 

in terms of the first-order simultaneous stochastic equations as 


^11 = ”1 

+ ^11^12 ^12^2/ 

^1+^2 " 'll 

yi2 = 

+ ^11^13 + ^11^11 

^12 

• • 

jm 

# ♦ 

> y 



= + ^11^1, rj^-2 

■*■ ^11*^1 r^ ^l*r^-l 

yi.r. 

= ni + 

+ ^l’'^?/r^+l+ ^1 rj. 


'^2»rj^+l = "^2 ^21^1r^ + ^22^2^ r^+2 ^ 2/ r^^+l 

^2/rj^+2 "" ’^2 ^22^2»r^+l ^22^2/rj,+3 ^ 2* r^+2 

• • • • • 

# • • • • 

^^"^1+^2"^ = ^2 ^2 2^2/r3L+r2-2 ^ ‘® 22^2# r^+r2"^ ^2»r^+r2-l 

y2.ri+r2 


= r)2 + P22^2.r^+r2-l ^21^11 ^2/r^+r2 



21 


These equations cam be written as 

^li ^ '^1 ^ ^lk^k/i-1 + ^ ^ii " ^ = 1#2. ' 

(2.2.2a) 

k = 2 where i = m = 2 when i = r^ * otherwise k=in=l 

and Yi, ^ = Yv ^ 

■^k^o k# rj^+r2 ' 


^2j ~ ^^2 P2k^k/j-l + ®2myrn»j+l + '^2j' ^ ‘ ' ^l'*‘^2 ' 

(2.2.2b) 


k = 1 when J = rj_+l. m = 1 when j = ^1 + ^2 ' otherwise 

k=fn=2 and y„ ^ _ , = Y^ i • 

■^m/ r+r2+l -^m^l 

In the above equations is the expected row yield of 

crop k in the absence of competition I and y 2 j are the 

random variables corresponding to the yields of ith and jth 
rows of first and second crop / respectively t is the 

competition coefficient between any two NN rows of crop k J 

/ k/&n= 1/2 represents the competition of crop m on crop k / 
'^li '^2j random errors and assumed to be independently 

distributed as N(0/Uj^) and N( 0 /CT 2 ^ / respectively. 

For convenience equations (2.2.1J can be expressed in 
matrix form as 


(i- 3 )y = 17 + u / 


( 2. 2.3) 


where v ' Y snd u are (rj^+r 2 ) component column vectors 
and B is ( square matrix given by 

2 = ’^l'^2 ^ 2 ^'' 1 = ^yil'yi2'*'*'yir^'y2/r3^+l"'y2/ 
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ion s rnsy oo wnitton 


( I-rij Y = PjLi + £ 


( 2 . 2 . 8 ) 


where Y = ( Y, / Y^; ' , 


2 ' ' 'i - = ( 




(rj_ i)rj_ 


^2 ^12 


P;i^ O 


H zz 2 


* P 


^1 ®21 


P]^ ~ ^k/^^l''‘^2^ * k = 1/2- 


(r2-l) ^22, 


^1 ~2 ^ Equations (2«2*7a/b) are exact and there is 

no need of approximations (2.2.6a, b). In practical situations 

and r 2 ^re likely to be small- 

Alternatively / the H matrix can be easily obtained directly 
from the parameters that quantify the total competition/ i.e. 


a'bi = 


2Cr^-l)8j,^ + 28^2 


2321 + 2 ( r 2 - 1)^22 


= H'aI / 


(2. 2.9) 


where 


= ( 1 / 1 . 


(lx( rj^+r2)) 


The equations for average yield per row are given by 

(l-2(rj-l)r;ip32)Y2 = ^ , 

where Y^ = Y^r^ and = t^r^. 


(2. 2.10a) 
(2.2.10b) 
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The total number of NN rows for r^ rows of first crop are 2r^» 
among which 2(rj^— 1) rows are of the first crop and two are 
of the second crop. 


In the case of MC the yield of ith row of crop k can be 
represented by 

^ki '^k ^ll^^k/i-1 ^k/i+l^ + ^i' ^ = l'2#..#r^+r2 * 

( 2 . 2 . 11 ) 


The corresponding equation for the total yield is 
’^ko = ^k + ^^kk^ko + ^ko ' 


^l'’'^2 

where Y, = S 
i=l 


^ki ' ^-ko 


^l'^^2 

2 

i=l 


ki 


The above equation can be written as 
“ ^^kk^^ko = ^k + ^ko' 


( 2 . 2 . 12 ) 


The extension to the plot with b repetitions of (r^+rj) 
rows is straightforward. The equation given by (2.2.8) holds 
good with the following changes ) 


N-, 


N, 


^k - ^ki ' '^^k = ^k\ • 


Prom equation (2-2.8) 

* E(Y) = (I-H)“^ PP , V(Y) = (I-H)"^ VP^/^ d-H')"^ 



(2- 2-13) 
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whe re P 


1/2 


: ° 1 


0 Vp, 


V = 


Vi 0 


‘ O V, 


and 


For the validity of equation (2.2.B)# (i-h) should be nonsingular# 
In the case of MC 

^">SV‘^ko> ' Vl-2?kk)‘^ • 

For example# in 212 row arrangement # i*e. r^ = r^ = 2# 
the first - order simultaneous stochastic equations # as 
described by (2. 2*1) # are as the following ~ 


yii 

= 

-h 

^ 11^12 

+ 

^ 12^22 

4 - 

^11 

^12 

= 


^ 11^11 

+ 

^ 12^21 

4. 

^12 

^21 

= ^2 

+ 

^ 21^12 

+ 

^ 22^22 

-f 

^^21 

^22 

= ^2 

+ 

^ 22^21 

+ 

^ 21^11 

-f 

^^22 


i* e. 

(l-3)y = T ) + u # 

where y = (Yh' yi2' Y 2 I' ’^22^ ' ' ‘ ' ; 


and 


B = 


0 

Bll O 
0 8 


^111 


21 


'12 


0 


O 


a 

H -7 O : 


P 2 I ' ° 


22 


0 



27 



Pll 

1 

^11 

^12 

ha": 


®11 

Pi 2' 

Here A' 3 = ! 

1 

1 

1 



i 





' Q 

^ 21 

L. 

fi 

^ 21 

^22 

022 i 


L 

^ 22 



—1 






i ^11 

^12 

' 



I 1 

0 0 1 

where H = 

^21 

^22 

i 

1 

A' = 

( 

D 0 

1 

1 1 


110 0 , 


10 O 1 1' 


The equation for the total yield is 
(I-H)Y = P 1L + £ , 


where Y^^ 




u 


ki 


The model (2- 2*8) also holds for certain situations viiere 
the interraw distances in MC 1 and MC 2 are different • For 
example / let ^22 interrow distances in MC 1 

and MC 2 # respectively * such that = ^2^22 ' where 

and ©2 are integers* The area occupied by rows of crop 1 
is equal to the area occupied by b 2 rows of crop 2- Let us 
consider an IC system where every set of b^^ rows of crop 1 
is replaced by a set of b 2 rows of crop 2 in MC 1 and vice 
versa • In IC the interrow distance between the crop 1 and 
crop 2 is ^ ^11+^22^ only consider the arrangements 

r-|_,* for which ^ 1 ^ 2 ^ ^2 integer* In this case the 

model C2»2«8) holds with p^ and p 2 as th«j following 

1 _ \ r / 1- i_ “1 


rj^/(ri + : 


= r2/(b2b“ r^ + r^) « (2.2.14) 


ha' . 
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2- 2* 2 Second-order simultaneous autoreqressive model 

Consider the situation in which yield of a row depends 
linearly on the yields of first and second NN rows* In 
• ^2 total number of neighbours for 

r^ + ^2 rows of crop 1 and crop 2 are 4(r2 + ^ 2 ^* 
Componentwise / the first and the second-oirder neighbours 
of the same and different crops are as following * 

Crop 1 Crop 2 



First-o rder 

second-o rde r 

First-o rder 

second-order 

Total 

Crop 1 

H 

1 

CM 

2 ( rj, - 2 ) 

2 

4 


Crop 2 

2 

4 

2(r3-l) 

2t r2-2) 

^^2 


The matrix A'b 1 can now be expressed as 


a'b 1 = ka' 1 = 

2^21 + ^^^21 ^^^ 2 “^'* ^11 ^^^ 2 ”^^ ^22 

where 

{2(r^-l)P2^^ + 2Crj^-2)p^^^ } {2p^2+^^12 ^^2 

H = 

- ^2^21 + {2(r^-l)p22+2^^2-2^^22^^ ""2^ 

is the competition coefficient between y^^ and 
k/m = 1/2* In the terms of total yield the equation is same 
as (2* 2.8) except for H which is replaced by the one given 
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ve. In some cases it may be realistic to assume some 
ationship between first and second-order competition parameters 
h as 

, p<2) ^ p(2) ^ 

;re a.^t and a^ are known constants- 

sn the total number of competition paraimeters to be estimated 
lains same as in the case of the first-order model • For example 
2*2 row arrangement 

O 
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B = 23 
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3 SOME ALTERNATIVE MODELS 
■ 3-1 Condi t io nal mo del 

Following the notation used in section 1-2-2 # let the 
>nditional distribution of be normal with conditional mean 

n 

E(Y. I all others) = + 2 c. . Y. # 


( 2 • 3 •!) 
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The joint distribution of Y is multivariate normal (3esag/1974) 
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The vector of yields y follows a multivariate normal 

distribution with 

E(y) = (1-3)"’^ V and V(y) = (l-3) ^ S * 


( 2. 3* 4) 


where 


diag (E f 


(ri+r^lxCr^+r^) 




= diag • 


(2.3-5) 


Using the methods similar to that or simultaneous model the 
distribution of total yir-lds of rows of crop 1 ana rows 
of crop 2 is a bivariate normal with 

eCy) = Pt and V(y) = (I-H; ^ PV . 

The necessary and sufficient condition for the validity 
of the above model is PV should be symmetric and 

positive definite- The symmetric property of the matrix 
v'^P^^CI-H) leads to the condition 


12 of =f>21 "2 


(2.3.6) 

This condiaor, restricts the use of the conditional model 
considerably- 

2-3-2 Moving average model 

The moving average formulation (1-2-10) in terms o^ the 
individual row yields ^ is 

y _ Tj + AU + u i (2-3-7) 


! ^11 ^12 


where 


. ■^21 ^22 
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In terms 

of the 

■cotal yields the model 

is 






Y : 

r; P M -f ( I-f-P) U # 



( 2. 3«8) 





2(ri-l)r-^ "ll 


“12 



where 

? 





t 






2(r,- 

^2 °^22 



Under the 

same 

assumptions about s as in simultaneous 

model 


Y follows 

a bivariate normal distribution with 




E(Y) 

= P u 

and 

V(y) = (I+F) 

V p^/2 (i+f') . 

(2.3-9) 



In experiments where competition plays a dominant role it is 
more realistic to consider the dependence of on the phenotypic 
values of the neighbouring rows rather than on residuals - Numerical 
iterative method is more complicated for estimation of the 
parameters in MA model in comparison to simultaneous model since 
inversion of Ci+F) is involved in each stage of iteration. 

2*3»3 Regression model 

One can also think of using the regression model by assuming 
that the yield of any row depends linearly on the true but 
unknown expected yield iVy) of the neighbouring rows in the absence 
of competition. Following the same notation used earlier the row 
yields in (r 2 ^*r 2 ) arrangement for the first-order dependence can 




where P*. and are the same as defined earlier 
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The equation (2-3>12) is similar to the linear mixture 
model (Schef fe' * 1958) wit±i additional factor 2ir-^+r^J 
This additional factor represents the degree of intimacy of 
the tvjo species* For example/ it is unity if the species 
are perfectly intimate / and this haopons when r^ = =1* 

The coefficients effects of proportions 

and P represents the effect of rnc; i-ow arrangement* Kelder (1963) 
suggested the use of mixture models for analysing the competition 
experiments such as random mixtures involving grass species* 

Prom equations (2-3*10a/b) the mean and variances of Y are 
E(Y) = CI-K)”^PM and V(y) = . (2*3-13) 

Here the expected value of Y is same as in the case of 
simultaneous and conditional models * but the variance differs 
and is free from competition parameters* Estimation of p 
coefficients by least squares in this case involves minimizing 
the quantity- 

(Y-(I-HJ“^ j ' P^/^ VP^/^(Y-(I-H;~^ PP) * (2*3*14) 

2* 4 DISCUSSION 

Let us consider a MC experimental plot with N plants- and 
suppose the yield- - of the ith plant can be represented 

as 

0). = 5 + ^ W . . (i) . + e . - i = i-2-***- N - (2*4*1) 

1 j^i J 

where 6 is the expected plant yield in the absence of 
competition- B* is the common competition coefficient and 
e^'s are errors with usual assumptions- V,^j is the weight 
assigned to the jth plant in relation tn i-hP n-t-h --.r.--. 



For 


which may depend on the distance between them- 

simplicity we take = 1 if ' ■^•ero otherwise# 

where d is some fixed distance . Let N plants be distributed 
over the plot in such a way that there are plants within' 
the neighbourhood d. By summing tne individual plant yields# 
the yield per plot can be expressed as 

r = N.s + (n^-1) p'y + £ # (2-4.2) 

N N 

where Y = 2 y , # t = S * 

1=1 ^ i=l ^ 

Since the neighbourhood area A is constant then for sufficiently 
large N # N^^-l ~ AP # where P is the plant density # i-e- -the 
number of plants per unit area- Hence from (2.4-2J we have 

E(Yi = N6/(l-gP) # (2-4-3) 

where P = p^A- In terms of expected plant yield 

E((u) = 6/(l-pP) = l/(a+bO) (2-4-4) 

where a = 6 ^ ' b = -P6 ^ * 

This is a well established as'yinptotic plant density 
relationship (Bleasdale and Nelder# 1960 ) Shinozaki and 
Kira# 1956) - Vandermeer (1984a# 1984b) discussed a similar 
t-ype of competition theory as above for developing the relation 
(2-4-4)- 

The competition model (2- 4-1) described above is the same 
as Meades (1967) model (l-2.ll)- Hence from this competition 
model we have obtdned well established plant density relationship- 



However# there is a basic difference /i. e. in competition 
model individual plant yields are to be considered for 
estimating the pa rameters whereas in plant density studies 
average plant yields or total plot yields are used. Moreover# 
there is a difference in the error structures which has 
implications in the estimation of the parameters- Proceeding 
with the other models such as conditional and regression * 
one does nor obtain the classical plant density relationship 
given by ( 2* 4* 4) • Thus# the use of simultaneous model has 
an advantage over the others. The simultaneous models (2«2.l) 
suggested in this study is a natural extension of the competition 
model ( 2 . 4 . 1 ) in MC- The procedure used to obtain the model 
( 2 . 2. 6) for total yield for IC is same as the procedure used 
to obtain the classical plant density relation (2.4*4). The 
model ( 1 . 2 . 2 ) suggested by Vandcrmeer (19h6) in terms of plant 
yields is similar to (2.2-1) with higher order competition 
coefficients bur it lacks appropriate stochastic formulation* 
Kempton (1982) and Besag and Kempton (1986) also used the 
simultaneous model (1.2-12) to incorporate the interplot 
competition in varietal trials- 


The restriction 
model is a serious 

P 21 

to assume 3p2/P'>l " 


drawback in IC 


in the case of conditional 
situations* In general 


of different degree and there is no reason 

The conditio niil model may be appropria 
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in MC situation where there is only one competition coefficient 
involved. There is hardly any difference between simultaneous 
ano conc^tional formxilations when competition coefficient 
is small t as has been pointed out by Besag and Kempton (1986) . 

In IC situations / however# both the models are different* 

MA model in spatial process is an alternative way of 
representing the simultaneous model vdth reparameterisation 
of the coefficient matrix, riowever# it is difficult to interpret 
the parameters in the MA model whereas in the simultaneous 
model they have a natural interpretation. We have obtained 
the regression model (2.3.11) by making the yield of any row 
depending upon T) # the expected rows yield in the absence of 
competition. We have seen that this is same as linear mixture 
model (Scheffe*/ 1958) with an additional parameter representing 
the degree of intimacy between the species* As suggested 
by Nelder (196 3) this type of model may be appropriate in 
analysing the e:cperimental data from random mixtures involving 
grass species* 

In the IC experiments the competition plays a dominating 
role* Thus# it is more realistic to think that yield of any 
row depends upon the actual yield of the neighbouring rows 
rather than on their expected values# as in the case of regression 
models* Hence simultaneous model appears to be more appropriate 
in IC situation* For using the mo'^el (?.2.8) it is necessary 
to estimate the parameters through an experiment of reasonable 
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size* These aspects are discussed in the next chapter- 
In the proposed model the yield of any row depends linearly 
on the neighbouring row yields* In case of grain yield 
this assumption may not be entirely valid* In such situations 
it may be more realistic to consider yield of any row depending 
linearly on other characteristics sucn cis vegetative growth 
of the neighbouring rows* However# the model considered 
here may not be unrealistic in case yields of the neighbouring 
rows are proportional to the vegetative growth. Normally 
p coefficients will be negative due to competition for resources* 
However# there may be situations where the interspecific 
competition coefficients ?‘ 2_2 ^ "^l positive# i. e» the 

presence of another crop is beneficial to the first crop* 

This has been more appropriately termed as cooperation* 

The model (2*2*8) can be used in arriving at the strategies 
of the row arrangements for achieving the maximum output per 
unit area. This is discussed in the fifth chapter in more 
detail* In practice# however# there may be ouher factors 
such as plant density and different fertilizer treatments 
which are of considerable Importance. The simultaneous model 
suggested in this chapter for analysina the IC experimental 
data involving row arrangements can oe extended without 
difficulty to these situations* For example# the model can 
be extended to study the effect of plant density on crop 
yields by choosing p coefficients as appropriate functions of 
intra and interrows distances* Some of these aspects are 
discussed ih the fourth chapte r- Introducing different genotypes 
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of the crops or more number of C3X)ps will increase the 
number of competition parameters* These aspects are not 
included in the present study. 



CHAPTER III 


FITTING OF SIMyLa>.NEOUS MODEL 
3*1 INTKDEUCTION 

In this chapter the details of estimation of the parameters 
in the first-order simultaneous model when the data are collected 
from a replicated field experiment / is discussed. Since the 
least squares estimators are inconsistent for simultaneous 
models' as pointed by Whittle (1954) and further discussed 
by Ord (1975) f maximum likelihood method (ML) is adopted in 
estimating the parameters- This method is discussed in section 3.2- 
Since the estimating equations are nonlinear in parameters 
iterative procedure for obtaining the estimates is developed. 

The ejqgressions for large sample variance-covariances of the 
estimators are derived in section 3*3. For eliminating the 
systematic fertility differences block parameters are introduced 
and the estimating equations are derived in section 3»4 • In 
section 3*5 the performance of the estimators are examined by 
some simulation studies* 

3.2 estimation OF THE PARAMETERS 

3.2*1 Maximum likelihood estimators 

Consider a two component IC experiment with t row arirangements 
and two MC plots arranged in a completely randomized design . Let 
there be N rows in each plot with b^ repetitions of r^j^ • 
rows of component crops (i = l#2/-.*t)* Each treatment is 



41 


replicated K times* In section 3*2 it has been shown r hat the 
plot yields of the component crops in IC form the following 
bivariate sets of simultaneous stochastic equations • 


(1 - 2 ( 


^li 


-1) r- 


-1 

ii 


‘ll^-lij 




^^2i®12^2ij 


+ e 


lij 


(3- 2-la) 


(1 - 2(r2^-l)r2j_ P22^^2ij ^ ^2i^2 + ^^li^21^1ij ''' ^ 2ij 


(3- 2.1b) 


1.2 


# • • • # 


j =1.2.* 


K 


N. 


where N 2 ^ = N . 

. ^kij ^kij yield and error associated 

with the plot in the Jth replication of the ith treatment of 
crop k 't is the error associated with the mth row in 

the jth replication of the ith row arrangement of crop k* 

Based on the assumptions about the errors, ns * as mentioned 
in section 2*2 . Y^j = (Y^ ij'W' follows a bivariate normal 
distribution with mean and variance as 


E(Yy) = (I-H^)"^ i V(Y^ j = (I-H^)“^ pJ^/^VP^/^(l-H^) 


( 3 . 2 . 2 ) 


where 


-1 


^^li"^^ ^li ^11 ■^2i ^12 


r;^ 3- 


1 


= 2 


^21 


^^2i“^^ ^2i ^22 


diag ^ P2x^ ' ^ ' ^k = ^^k ' 


Pki “ ^ki/^ ^li'^^2i^ ' ^ “ 


1 . 2 * 
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In the case of MC 

' ''(^koJ> = Vl-2Skk>'^ • (3-2.3) 

where is the yield of the plot frora the jth replication 

in MC of kth crop- 


-21og L = constant - 2 K[s log iG^i + log { (1-2^ 1-28 } ] 


22 ' 


+ K(t+1) logCV^-V^; + S S(G^ - m)' V^Ig^ Y^ - jl) 


1 3 


t E [ v : 


-1 


{ (1-20, ,) Y 


kk ko j 


ij 

- M , - 


-iJ 

(3*2.4) 


where G. 




V. 

i 


vpTI , 

1 


By assuming 8^^» 0^2 ' 3 22 ^21 ^ estimates 

of }J-^f U- 2 ' snd V 2 are obtained by equating partial derivatives 
of -21og L to zero* and they are 


K(1 I = E E I ) 


* < l - 28 ll ) E Yjoj < 


( 3 . 2*5a} 


1 ' 21 - "2 i j 


K(1 -f E p„j H, = E S I (l-2Cr2i-l)rJ^ 922^^13 * 3ji Yiij ! 


+ (l-2Bjj) E 

- 1 . 


(3* 2- 5b) 

- 1 . „ A ,2 


K(l-t.t)V^ = E S J ^^"2(rj^j_-l) r^^ 83 _^)Y^^j - 2r2i8i 2^2ij"^lih ^ 


A ,2 


+ S I ( l - 28 ii ) Y^^j - 


( 3 ♦ 2 • 6 e ) 


K(l+t)V2 = E E p^i { (l-2(r2i-l) 322^^2ij ~ ^^li^21^1ij “P2i'^2 ^ 

1 . J 


+ E { (1-2R^^) Y^qj - txo ^ ^ 


(3.2.6b) 
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By substi taring these values {3* 2- 4) reduces / in terms 
of competition coefficients * as 

f(B) = constant - 2K [ S log IG^i + log f (1-2B3_) (1-2^22^ ^ ^ 

A A 

+ K(t+1) log (V^.V 2 ) » (3*2-7) 

where p' = (B^l ' ^12 ' ^2'^ * ^21'** 

The ML estimate of 3 is the value corresponding to the minimum 
value of rhe function (3*2*7J. 


3* 2- 2 A Nimerical Method for finding ML Estimates 

From equation (3* 2*7) the ML estimate of 3 is the value 
corresponding to the minimum value of the function 


f^(e) > log MS E p-1 -2(t+l)-^ [ Z log 

K 1 j J 1 

+ log U-2i3^^; (l-23,p) ] (3*2.8) 


22 ' 


i = 1» 2 


/ ♦ • • # U # 


j = 1» 2 


f • * • f 


K ,/ k =1.2- 


where fj^CB) = (f(B) - constant) /K(t+1) y 

®kij ” -2(r^^-l) ^^kij ~ ^^mi ^km ^mij ” k . 


( 3 . 2*9) 


m = 2 -vmen k = 1 m = 1 when k = 2 # 


®koJ = '^koj - *"h ' (3.2.10) 

“ Z * ^kl “^kl ^kk'" ^ ^12 ®21 

iC 

( 3 - 2 . 11 ) 


and is as in (3.2.5a.bJ. 


•p> 
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The coefficients are obtained iteratively* The estimate 

at the (n+l)th stage of iteration is 

(n+l) oCn) .,-1 


where ^ 4^ 


f^Cp) = ( 3 f 2 _/ 9 0p) / H is the Hessian matrix i*e. H = ( 3^f^/3Pp 3 P^) » 

p# q =1/2# 3' 4 • 


The elements of fnCS) and H are given below- 


9f^/3 8 


= 2 [ E S + S C e^oj^^koj^p j ®3cij + ^ ®koj 


2 3-1 


*"1 t T~* / ^ \ / ^"1 \ ^1 


-2(t+i; { S (d.)' d7 - 2(l-2p ) } 

a. P 1 q 


k = 1 when p=l/2 / k = 2 when p = 3/4# and P^ = P^^ when p = 1 J 

q 


P = P -5 when p = 3 / otherwise zero 




^kij + 2K- p^^Cl+I p^^)- . For p = 1- k = 1 

, ) 

.[ E E I (r^i-l)r-t ) + E ] For p . 3. k = 2 


1 J 


L-2^ml + 2K-Xi<l + \ Pkl' "" ' ^ E r^J ) For P-S'l'-l 


T.-l 


-1 


For p=4/k=2 an 
m=l 


=9e.^,,-!/ 3Pr-> can be obtainea from Ce^^jip by replacing 


'koj p 


'ko j 
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;he first term by / when p - 1/3 and by zero when 

• = 2/4* 




-2(Pj,^P2ii"^(r^^-l)r"^ U-2( r'^0pl ForpJ..k=l and in.2 

For p=3/k=:2 and m=l 




For p=2/ q=4 
For p=4/ q=2 


rhen p/ q = 1# 2 


Kp.q) = S 2 + I [| s + ®loJ>p>- 

{2 E + ^ ^®loj^ ^®loj^q 

+ 2(t+i;“^ E { (d^)' (d^J ' / } + 4 K gCg) / 


«?here g(p) = (l--262_) when p = q = 1 ^ and zero otherwise • 
?or p =s q =r 3^4^ ^lij ^loj replaced by ^2 ±j 

^2oj ^ respectively/ in the above expression to get H(p#q) # 

iere gC3) = (1— 20^) when p = q =t 3 and zero otherwise- 

•^hen p=:l#2 / q = 3/ 4^ 


H(p,q) = 2(t+i)-^ [ E (d^);(dy^ / 4 + 1 / \] 


■There 
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r ^ 

ki for p = 1. q = 3 

K 

^2i^ ^ for p = 2/ q = 4 

O otherwise. 

S. 

smaining elements of H can be obtained from the property of 
Tnmetry. 


3 VARIANCE-COVARIANCES OF THE ML ESTIMATORS 

The expression for the large sample variance-covariance 
itrix of the estimators f discussed in section 3*2 # is derived 
►re* The log-likelihood is 


I s log L s constant + K [ S log d^^ + log [ (1-20^2^) (1-2022^ ^ ] 


- 5 i °9 ‘ A - Vj ) - Is E ( G ^ - M )' v - Rg ^ y ^^. p 


i J 


-ij 


- 3 ^ [ ''k^ ^ ‘ 'i-2ekk> ’^koj - “k > ^ 


lere d^^ is given by (3 •2*11) 

I the case of IC let 
,-l 


E ( Yj ^ j ) = g ”-^ li = 


-1 


^i ' 


^ = -i -i + ' say / i = 1. 2# • . ^t # j = 1# 2- 


-ij -ij 

®i = ^® kmi ^ ' ^i 
(2x2) (2x2) 


^^kmi^ ' -i ~ ^^li'^2i^ ' 


1 the case of MC let 
koj' = ^k'-"""^kk^ ' ""ko ' '''"koj 
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Simplified expressions* after differentiation and taking expectations 
are as the following • 
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elements of the symmetric matrix B are as the following 
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The ez^ected values B(l* 4 ) and 3(2/3) are given by 


:(B(p/q) = _K [ s r (d^)' (dj_)' d”^ ] ]. 


For p =t 2/4 

E(B(p,p)) . -K[ s Udi); -4V^- E p,^^ 


Forp = 2/ ]c=l/ and m = 2 .* forp = 4/ ]c = 2 and m = 1< 


3»4 FITTING THE MODEL WITH BLOCK EFFECTS 

In agricultural field e>qperiments blocking will usually 
eliminate large scale fertility differences* By considering 
the arrangement of the rows in randomized blocks/ the 
simultaneous equations after introducing the block parameters 
are 

"^lij = Pli^^l+^l ^^^li“^^^li ^11 nij ^^2i ^12 ^2ij ^lij 

( 3 ♦ Ih) 

^2ij = P2i^^2+“2j^ 2(r2i-l)r2^ ^22 ^2ij ^^li ^21 ^lij ■*" ^2ij 

( 3 • 4* lb) 

where ct^ j (l- 2 P 2 ^j_) and “2j^^''^^2 2^ effects of 

the jth block on the first and second crops* The expressions 

that estimates b- 2 ^ and i ^2 (3*2*5a/b) will estimate and. 

/ where a = E a, ./K / k = 1/2 . The estimating equations 
2 2 ^ J 

for error variances and V 2 are 





tr 

< > 0 





ijj 

d 

(D 

M 

hh 

d 


H* 


• 





(D 




K> 

CA 

P 

P 

(t 



■f 


• 

CD 

0) 

d 

d' 

Hi 




H 

0 

(D 

pi 

CD 

P 


P* tr-j 


o 

<t 

pi 



d 





H- 


<> 

0 

Q 





0 

H- 


0 

ft 


to 


T:> 

d 

ft 


M 

H- 


H* 




Q 

CD 


0 





Ia) 


h 

CD 

d 


.<~V 


H- 

• 

P 

CD 

ct 



p 


cn 

K) 

(t 


H* 

d 


lu 



• 

H- 


ft 

0 


H* / 


ft 


<1 

0 

p. 



P 


0 


CD 


0 

tr 


M 



% 


O' 

d 

(D 


C-#, 


cr 



CD 




v.^ 


CD 

h 



Q 

3 





(D 


0 

0 

p. 


II 


h 

3 

0 

tr 

CD 

d 

1 


1 

CD 

P 


ft 

Hi 

R* 

fvj 

p. M 

to 


H' 


CD 

Hi 

3 

. H . 


d, 

H 

P 

0 

H* 

H- 

H- 

H* 1 


to 1 

OJ 

K 

D' 

d 

O 

^3 

p. H 


p. H 

0 

Cft 

rt 

CD 

H* 

CD 


H 


CD 

6" 

CD 

Pi 

CD 

Pi 

"CD 

r 

TO 

p. 

3 

H* 


d 



fO 

P* 


u) 

d 


ft 

hh 



to 

tr 


H* 


(A 

O 


H 


KJ 

CD 

d 



H 


to 

K 


X 

lO 


p. 


H* 

H‘ 

to 


o 


tjJ 

(A 

0 

p. 

1 

p. 

T 

c^ 

rt 

• 


0^ 

L-J. 

pj 

Cj. 

^ V 


d' 


(A 

ft 





ft 

CD 

• 

CD 

CD 

vw 

H , 


p ^ 




3 

p. 

4- 

M 1 

4“ 


i tl 

3 

pi 

CD 

d 


p. Pi 


L>j, 

0 

p 




pj 


— ' 


h 


tr 

P> 

d 

1 

TO 

M 

0> 


CD 


(A 

vQ 

fv) 

fO 

1 

w 

t; 

</i 




XD 

i\^ 

NJ 



(t 

ft 


ft 

K) 

w 

CD 

iO 


tt 



d 



p* 

H* 


B 


• 

CD 


K 

H 

< 


a> 

h3 




to 

“ — 


H 

ft 

tr 

« 

S 

K 

P- • 


P 

D 

(D 

(D 




f~i* 

H* 



W 




0 


0 

H* 



(D 


(D 

Ui* 


Ui. 


O) 

% 


CD 

CA 

• 




• 


p 


ft 



% 


K> 

pi 

p 

O 

p. 




Oi 

• 

H- 

<t 

CD 

3 




• 

VO 

0) 

H- 


0) 




»rv 


0 

0 

ft 

(t 



u» 

• 


P 

d 


(D 

* 


• 


CD 

W 



CA 




cw 

d 

(A 

ft 



« 


« 


Pi 

CD 

0 



U) 


OJ 

tr 


Pi 




tr 


CD 

• 








«w» 




1 






H D- 






CD 


ft 



ft 

•f 

1 


-t 



4* 

rt) 





M 

H* M 







V 



< > 



<> 

H 



M 



p 











li 



!l 


•f 

1 

p. M 

+ 

t 

p. iM 



N) 


L_». M 

t\> 



L^. M 


Ui. M 

h 

O. M 

-f > 

p 


-1 

li 

"d 


-1 

2i 

"d 

p. 



^v) 1 



h- I 


H 

TD 

p, pi 

p 

xo 

1 

i 


f 

M 


1 




fO 

H 


(U 

fO 


II 

XD 


tr^ 

TD 



M 

K 

H 

H 

k; 


M 

ro 


1 

p 

CO 

H 



p. 


■V^ 

p 

I 



Lj. 



iwl. 

ro 




hj 

p 


.<~s 

H 

0 

1 

f\> 

0 

1 

h 

1 

i-J. 


p. 

i-i. 


H 

fv) 


tJ 

1 


n3 

p. 


1 

fO 

p. 

1 

h-* 

1 

ht 


H* 



H* 

H 







N— ^ 

p. 

p 

■p 

fU \ 

T: 

■p 


1 

LSJ ^ 

^ N 

. H* H 


h-* 

d 

H 

+ > 

4- ^ 


+ / 


p , 


R 

R 

T» 

P 

-f ^ 

p. p 

li 

K) 

N> 

H 




H 1 

Lj, 


H 

Lj. 

P 

XD 

h»* M 




V-^ 

p 


a> 




«.w 

t-j. 

p 

M 

(.W 

s-»»-» 


to 



M 

N) 

CO 

N) 


'-w 


'w-' 



1-^ 


K> 

K 

>-< 

% 


Ui. 

% 


p 

M 






H- 

H* 






Uj. 

Uj. 

<r-^ 







CO 



(Jj 







* 











• 



• 




N) 



to 












11255S 


3*5 MONTE carlo SIMULATION STUDY 


A simulation study is carried out to study the performance 
of the ML estimates of the parameters for an experiment of the 
moderate size- The Monte Carlo experiment consists of 36 plots 
arranged in three randomized blocks- Among 12 plots in each 
block 10 plots are with row arrangements of the first and the 
second component crops as 1*1- 1*2- 1-3' 2^1- 2^2- 2-3/ 3*1# 3*2- 
4^4 and the remaining two plots are of MC. Values of the 
parameters in MC are taken as ^ ' ^2 “ ^11 “ ^22 ~ -0*4- 

The values chosen of the parameters are such that the expected 
yields are in the range of vAiat we find in sorghum and pearlmillet 
MC« Three situations of interspecific competition are considered 
by choosing ^21 

i) (-0-1# -0-1)# (ii) (-0.1# -0.4;# and Ciii) {-0.1#-0.6)# 
respectively. 


Data for 500 experiments are simulated with and Vj such 
that the coefficient of variation (CV) are 10% and 20% for 
two MC* s - respectively# which is the usual range to be observed 
in the field experiments* The number 500 is chosen as the 
changes in the average values of the estimates are considerably 
smaller for higher numbers* The IMSL subroutine (GGNRM) 
package is used to generate the pseudo bivariate normal variates* 
Block effects are assumed to be zero- The ML estimates of ^*3 
have been obtained by solving the nonlinear equations - using 
Newton - Raphson iterative method described in section 3*2# by 
including the block parameters* The starting values are obtained 
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by simplex subroutine (Nelder and Mead,1965J of NAG. The 

results for three different cases considered are given in 

Table 3*1 • The ML estimates obtained by the Newton -Raphson 

method are close to the starting values given by the simplex 

subroutine. The convergence of iterative process did not 

create any problem except being slow in about one percent of 

the cases- Convergence is said to be achieved when the 

absolute values of the partial derivatives of the p coefficients 

■“3 

are less than 10 and the successive values of fj^(p) had a 
relative change of les^ than 10 . The average number of 

iterations required to converge is found to be approximately 
three. 


Relative bias and MSE for the estimates are calculated as 

A 


bias CB^) = . 

A A 

relative bias x bias (6 


(3.5.1) 
(3.5. 2 ) 


R 


where ^ 3 ^. = / R , R = 500 

‘Pkml - 

1=1 


(3.5.3) 


The values of bias and '^MSE for the three interspecific 
competition situations at CV 10 ^ and 2096 are given in 
Table 3*1 . From the results (Table 3*1) it is evident that 
the bias and the '^MSE are dependent on CV and also on the 
' values of the parameters to some extent - This is ejqpected 
from the e 2 Q 5 ressions of lar^e sample variance—" variance matrix 
in section 3-3. Bias and '^MSE are considerably large at 2096 



Table 3-1 Results of simulation study for three sets of parameters values at 10 94 and 


5 3 



cu o 
cn as: 
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than that at 10?6 , as expected* The parameters and are 
estimated with a relative bias of not more than 3 % .in the 
case of Competition parameters the relative bias varied 
upto 15% when CV is 10% and to 27 % x,hen CV is 20 % . In 
most of the cases bias is less than onetenth of '^MSE- In all 
the cases the bias contributed only a small proportion of 
the MSE* For most or the cases the variance of the estimates 
is at least 99% or the MSE. In general it appears from this 
simulation study that the parameters can be estimated satisfactory 
at 10% CV except when the competition coefficients are very 
small. Vfhen the competition coefficients are small then the 
relative mean square is large • as expected. However<^ in such 
cases relative mean square error is not a very good measure 
of the variability* 

The effect of increasing the number of replications to six 
and nine has been examined* The results are given in Table 3*2 
for one set of parameter values at 10% and 20% CV* The bias 
decreases when the number of replications are increased and 
the variances decrease proportionately with the increase in the 
number of replications* This reduction is expected* The maximum 
relative bias reduced from 14% to about 10% and 9% at 10% CV # 
and from 21% to about 16% and 10% at 20 % CV when the replications 
are increased from three to six and nine ^ respectively* This 
suggests that when variability is large and competition coefficients 
are small* large number of replications are needed to estimate 


the parameters precisely* 



GV 


10 % 


20 % 


^11 = ^22 

= - 0 . 40 / 

^12 “ 

P 21 * 



^11 

^12 

^22 




Replications = 

Bias 

2 

X 10 

1 

- 1.51 

- 1. 41 

- 0.91 

'^MSE 

X lO 

1. 60 

1.55 

1.70 




Replications - 

Bias 

X 10 ^ 

- O . 20 

0.16 

- 1*33 

'^MSE 

X 10 ^ 

1.12 

1 . 0 ^^ 

1. 25 




Replications = 

Bias 

X 10 ^ 

- 1.17 

- 0. 39 

- u . 81 

'^MSE 

1/^1 

X 10 

0.96 

0.91 

0.96 




Replications == 

Bias 

X 10 ^ 

- 2.54 

- 2.00 

- 2 . ^0 

'^MSE 

X 10 ^ 

2 . 25 

2.19 

2 . 30 




Replications = 

Bias 

X 10 ^ 

- 1.10 

- 0.07 

- 2.45 

'^MSE 

X 10 ^ 

1.50 

1 . 47 

1.58 




Replications = 

Bias 

X 10 ^ 

- 1.69 

- 1.00 

- 1.82 


X 10 

1 . 23 

1.18 

1.26 



50 / CV = 10 96 and 20 % 


iLLj^ - 60 / 4 ^ = 


^21 


^2 

- 0.69 

99.79 

60. 69 

1.06 

112. 87 

98.13 

6 



- 0.95 

0.34 

83. 42 

0.77 

73.05 

31.13 

9 



. 0 ) . 40 

70. 68 

42.71 

0.58 

66 . 43 

54.13 

- 2.10 

161. 28 

162.07 

1.53 

155.78 

136.00 

6 



- 1.62 

48*45 

152.41 

1.06 

103.54 

95 .81 

9 



- 0.08 

94.67 

94.89 

0.79 

84*36 

71. 40 
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There is a close agreement between the large sample 
va riance— covariance matrix obtained through the expressions# 
derived in section 3*3 # and empirical variance— covariances 
obtained tn rough simulation study* For example # the variance- 
covariance matrix from simulation study ^Emp.V(^) # and the 

. . A 

asymptotic variance-covariance matrix ^ Asy. V(p ), for the 


competition parameters are given 

below 

for the parameter values 

^11 = ®22 "" -0*40 # 0^2 = ^21 " 

— 0 • 10 # 

= 60 < 

.00 # ^2 = 50.00 

and CV = 10 % . 


-t 


2-57 2-28 

-0.90 

-0. 61 


Emp, V(0) = 10~^ 2- 39 

- 0.97 

-0. 65 

0 . 5 . 4 ) 

— 


1.67 

1.13 


2- 46 2-11 

-0.78 

-O. 52 I 


A -2 2-19 

-0.87 

-0.56 

( 3 . 5 . 5 ) 

\ 

{ 

Asy. v(p) = 10 

2. 68 

1.49 ^ 

where 0 = 2 ' ^ 22 ' ^ 21^ 


0.99 



A similar type of agreement is observed for the other cases - 
3 . 6 DISCUSSION 

In the absence of exact tests# we shall use tests based on 
the asymptotic variances and likelihood ratio # which are only 
approximate for finite sample sizes- The ML estimators obtained 
here are consistent but this gives no idea of the magnitude 



of bias which are likely to occur for any given size of the 
experiment. The simulation study carried out here gives some 
idea about the magnitude of bias and approximations involved 

in using large sample expressions for variance-covariance 
matrix. 

In general there exists competition among the crops^ 
consequently tne competition coefficients are taken negative 
in simulation study- These Goeffici«nts are selected in such a way 
that the expected yields are realistic Willey (1979) describes 
three different situations in IC. He namef them as i) ^mutual 
cooperation ir the intraspecific competition of both the 
species is greater than their interspecific competition# 
ii) 'mutual inhibition' if the intraspecific competition of 
both the species is less than their interspecific competition# 
and iii) 'mutual compensation' if the intraspecific competition 
of one or the species is greater than its interspecific competition# 
and for the other species the former is less than the latter* 

The first type of competition is not unusual. The second type 
is rare and in this case IC is clearly not advantageous over 
M.C. The third type of situation is most common* Among the three 
cases selected for simulation the first one represents mutual 
cooperation# and the second and third represent mutual 
compensation - 

The choice of the initial values for starting the iterative 
procedure # discussed in section 3*2.2 # is important otherwise 
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there will be a problem in the convergence* The simplex 
procedure adopted gave satisfacto r** values of the estimates* 
the Newton-Raphson method further im.j roved these estimates 
to satisfy the required convergence criterion* The possibility 
of convergence to a local maxima has been examined in detail 
for some cases by giving different initial values# but the 
function always converged to the same value- This rules out 
the possibility of convergence to a local maxima* 

It is evident from the values of MSE that more number 
of replications are required -v^en variability is large or the 
competition parameters are small in magnitude* One can search 
for optimum designs such that the parameters can be estimated 
precisely through an experiment of reasonable size . 

If the individual row yields are available it would be 
interesting to examine the gain in efficiency in estimating 
the parameters* However# this will involve a lot more work 
in data collection and computation. In IC experiment it is 
a normal practice to record plot vields ana consequently this 


has not been pursued* 



CHAPTER IV 


SPATI^Ij model with varying population densities 
4 . 1 introduction 

Plant; population of the component crops is another major 
factor that atfects the total competition pressure per unit 
area / as has been pointed out in the first chapter- The 
knowledge about the nature of response to varying plant 
densities will help in arriving at the optimum density levels 
for achieving maximum yield- The relationship of yield and 
plant populations in IC needs a detailed investigation as 
in the case of MC (willey/ 1979j . In MC Willey and Heath (1969) 
reviewed the practical aspects of the relationship between 
plant density (plant population per unit area) and crop yield- 
The inverse power response curve suagested by Bleasdale and 
Nelder (I960) is the one which is widely used. Mead (1970) 
examined the theoretical properties of this relationship and 
also discussed estimation procedures* Gillies and Rotkowsky 
(1978) found that the maximum likelihood estimates are biased 
and to overcome this partially they suggested a reparameteri- 
zation- The plant density relationship based on the principles 
of competition in MC is discussed by Watkinson (1980). This 
is again a reparametarized form of the Bleasdale and Nelder 
equation- The parameters in this equation describe the 
biological process of the competition and hence ensure wider 
applicability- The equation developed in this chapter for IC 
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out to be a natural extension of this relation* 

The model suggested in the second chapter has been 
extended to include the effect of plant population at a 
constant row width. The estimation procedure for this 
extended model is discussed for field experiments in 
great detail. A particular case of the model in terms of 
plant densities only is also considered. 

4.2 DEVELOPMENT OF MODEL 

In section 2.4 we have seen in the case of xMC the 
classical yield density relation 

e(y; = N6 + pp eCy; , (4.2.1) 

is derived by giving equal weightage to the neighbours 
which lie within a fixed distance. Now ve examine the 
effect by varying interplant distances* As the distance 
between the neighbouring plants increases the competition 
decreases. Let us assume that the relationship between 
competition coefficient p and the interplant distance x 
be given by 

6(x) = Pq » (4.2.2) 

where 3^ and c are constants. There is an empirical 
evidence for this relationship (Vandermeer/ 1986) . Let 
the interrow distance / x^ # be constant # then 

-1 -1 
p_(x^x>l /i.e.Pccx 9 

and (x x) is the area occupied by a single plant. 


(4.2.3) 
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From ( 4 . 2 . 1 ) / ( 4 . 2 . 2 ) and (4.2.3) 

c, 

E(y) = N6 + P E(Y) , 

where = c+1 . The expected value of the yield per plant 
a)( = Y/N) / is given by 

E(m) = 6(1 - P^j'^ = U + b , ( 4 . 2 . 4 ) 

where a = 5 ^ and b = /6 . 

The general form or the yield -• density relation suggested by 
Bleasdale and Nelder is 

Ed/oc^) = a + pp*^. ( 4 . 2 . 5 ) 

The ratio 0 / b determines the shape of the curve . It is 
asymptotic if B / ^ ^ 1 # otherwise there exists a maximum at a 
finite value of P • Based on this fact Bleasdale (1966) suggested 
ff ss If which belongs to the family of curves given by (4.2.5) » 

Apart from the approximation E(l/co) ^ 1/E(ci)) * the equation ( 4 . 2 . 4 ) 
is a particular case of (4.2.5) when 0= 1. Approximation may be 
justified by ignoring the higher order terms and noting that o) takes 
only positive values. Hence the model developed above# based on the 
principles of plant competition in MC / belongs to the widely accepted 
Bleasdale and Nelder equation. The equation developed in this 
chapter ror IC is a natural extension of (4.2.4)* 

Development of model for 2.2 row arrangement 

We shall start with the competition model at the individual 
plant level in IC. From this^a model for the plot yields # 
including the plant density effects # will be developed. 

For simplicity consider 2 *. 2 row arrangement . The assumptions 
regarding inter and intrarow spacings in IC and AC are the 
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same as discussed in section 3 .?. Let y, , and v, ^ be 
the yields or rows 1 and 2 of the ~irst species • Y'23_ 
and Y22 yields of the two roves for the second 

species# and yields or ith and jth 

plants in kth and mth rows of first and second species # 
respectively # (k#m = 1 » 2 ) • We assume here that the 
yield of any plant within a row is affected by its 
immediate neighbouring plants apart from the iiranediate 
neighbouring rows- Let there be n^ plants in a row of the 
kth species* If the competition coefficient due to the 


neighbouring row is # then on the average each plant 

ejqjeriances a competition of to the neighbouring 


row- 


Let and 0)2ij 


be the yields 


an j th plant 


of the first rows of crop 1 and 2 » respectively# then# 


^lli = + Pl^^ll#i-1 + “ll#i+l' ^11^12/’^! ®11; 


( 4« 2 * 6h) 


'2lj = ^2 + ^2^“l2#j-l + “l2#j+l^ ^ 32 '^ 22''^’^2 ^21^12/^2 ®2lj 


1,2# n^#' j = 1# 2# . * • #n2 ' 


( 4 • 2 • 61d) 


where 5 ^ and 62 expected plcint yicads in the 

absence of intra and interrow competitions * * ^2lj 

are the errors assumed to be distributed inaependently as 
N( 0 #e|) and N(0#e2) ^ respectively* 
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By adding the equations corresponding to all 


a 


Hi 


and f»^21j ' have 


Yll _ + P12 ^22 "*■ ^11 ^ ^ 4- 2.7a) 

^21 = “1^2 + 2^2^21 + ^22^22 + ^21^12 + ^ 4- 2.7b) 


n. 


n. 


where 


>'ki = “W1 ' “h = ®iai • 


The a.bove equations can be written a< 


yil = 'll + Pii yi2 + 612 722 + “11 , (■S.2.8a) 

^21 = llj ■*■ ^22 ^22 ■*■ ^21 ^12 ■*■ "21 ' (4-2. 8b) 
where = hj. ' "ki - 


These equations agree with the equations formulated earlier 
for row arrangements (2.2.1). Now let us examine the effect 
of varying number of plants within a row by keeping the 
interrow distances fixed- From the average per 

plant yield of the first and second crop can be expressed 
as 


0 ) 


-1 


11 


= 5 i + 0)^2 + ^12 ^^2 ^'1 “22 ’^ ll ' 


(4*2*^ 9a) 


“21 ~ ^2 + ^^2 “21 ^22 “22 ^21 “^1 ^2 “12 ^21 ' 


,21 = 02 + ^P2 '"21 


(4.2.9b) 


where 


“km = W’" 


k ' 


^kl “ '^kl/”} 
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For a fixed row length the number of plants within 
a row is proportional to the average distance between the 
plants. Let x^ and x^ be the average plant to plant distances 
within rows Ox. the first and second species# respectively # 
and Xq be the fixed interrow distance. Prom equations (4.2.2) 
and (4.2.9a/b/ # we have# 


Cl) 


11 


= 6i 


lo 


_ 

* 0) 


11 


+ p 


110 


X -, 


'11 - 
(d 


12 


+ n^ n 


-1 


0120 


'12 ^ 




22 " 1 


(4.2.10a) 


-c 


Cl) 


21 


6 o 20 


20 


X. 


2 - 

Cx ) 


21 


8l 


-c 


20 


22 

G ) 


22 




n. 


-1 


-c 


0 21 ^ 




(4.2.10b) 


where c^^ # c^^ # k#m = 1#2 are some constants* 
Since the interrow disrances are fixed 


x~lo<t(x^ x^) 1 = « 

where plant density of the kth species- 

The above equations can be written in terms of plant 


densities as 


'11 


0) 


= 6^ + 23^0^1 “ll 0110 ^1 “l2 "2 ’^l 0120 ^2 “ 22 '’' '^ll ^ 


'12 


11 


( 4 . 2.11a) 


-2 '22 , °21 


= 62 + 23 ' P2 ^21 ”^ 0 220 ^2 ^^22 ^1 ^2 02 


(0 


21 ~ ^2 


20 "2 '"21 ^ ^220 


210 *^1 “12''' '^21' 


(4.2.11b) 
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The equations in terms of row yields are 

^11 ~ ^^10^1 ^11 ^ ^110 ^12 ■'* ^120 ^2 ^22 ^11 * 

( 4 • 2 • X 23.) 




721 = n 2'^2 + 2 B '20 ^ ^21 + ^220 '>2 ^22 + ^210 ^12 + 4 

(4- 2.12b) 

Similarly we can write the equations for the other two rows 


c c c 

Y^2 = "iSl + Y^2 + p/2 

( 4 ♦ 2 • 1 33 j 


^ O ^00 

^ ^ . n* t\ 2. ^ 


^22 “ ^2^2 ^^20 ^2 ^22 ®220 ^2 ^21 ■*■ ^210 ^1 ^11 ^22 * 


(4. 2.13b) 


From the above . equations in -cerms of total yields are 


c c 

= 2 n^ 63 , + Pj,^ + B'lO Pj 


^ 8^20 ^2 "l 


(4* 2*1 43 ) 


Y2 = 202^2 + 2B ' 2o P 2 ^ + 822 ^ P 2 ^^ ^2 + ^210 ^ 1 ^^ ^1 + ^2 ' 


( 4 . 2.14b) 


where + 7^2 ' " kl + '^k2 


Model for monocroppinq 


In MC the jth plant yield in the ith row of the kth species ^ 


“koij ' expressed as 
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in IC / we have the rtradel in terms of plant density as 
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The above equation is same as (4* 2.4). For simplicity we 
reparameterize and as 

4 n^ 63^ = 4n^ 6^ = n^^ P^^^ , 

where <^ 3 ^^ = 4 n 3 ^^ S 3 , ^ = p 3 ^ * P^^ is plant density 

when there are plants per unit length of crop k. 


Hence 

(4.2-18) 

can be 

written 

as 

Y 

■ko 

li 

0 

+ ^^ko 


+ 2Skko "'to ^ "to ' (4.2. 20) 

where 

^ko ■ 

8' 

^k ko 

' ®kko 

®kko ("ko ■ 


Generalization to rj^ '• r^ arrangement in IC 


From (2.2. 8 ) and (4*2.l4a/b) the equations for the total yield 
of r^ rows of ti 


■ 2 = ^2 “2 °2 


first 

; crop and 

^2 

rows < 


'"l 

P/ + 

2( 

-IJ r^^' 

2r;^ 

^120 2 

^2 

+ £3^ 


'^2 

^2 ^2 ^ 

2(r, 

,-l)r; 


8' P""^^ 

P 2 IO ^1 


+ ‘'2 


'11 


(4.2.21a) 


'22 


(4.2.21b) 


For the case of MC we have 


^ Iclc 

^k ■*■ ^^ko ^k ^ko “^kko *^k ^ko ■*■ ^kO 



68 


! can write this equation as (4.2.20) , by taking 


(r, +r^) n, • 

1 2 ko k 


e equations (4* 2* 21a/ b) can be written 


=■ Pi h % + 233^0 + 2 ( 13 ^-!) 


+ 2r"^ 0 d ' (4.2.22a) 


^2 ~ ^2 ^2 *^20 + 26 d 


20 ^2 ^2 ^^^ 2 “^^ ^2 ^220 '^2 ^^ ^2 


^210 *^1^ ^1 + ^2 ' 


( 4 . 2- 22b) 


Pkmo^ Pk = ^k / ^^1 + ^2^ 


juations (4*2*22a/b) can be written as 


o ^ 

-2(r3^-l) r'^ PjLio !^1 = Pi <il ‘'10 + ^^'2 h20 '*2^^ b+ h' 


(4. 2.23a) 


11-2^20 + 2(r2-lir;^ ^220 ’b ” P2 '^2 "20 + ^210 ^l+=2 


( 4 . 2.23b) 


urther this can be e 3 <pressed in matrix notation as 


(I-H)Y = P D + e / 


( 4 . 2 . 24 ) 


here Y = ' fo = * D = diag (< 3 ^'^ 
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Oslo'S!" - 2(ri-l)r-^ 


(I-H) = 


^ _'>-r~^ S d 

“-2 ^120 ^2 


-2r-^ 3 210 ^1 


21 


1-23^0^2^ -2(r2-lir2^ 3 330 ^2^^ 


In MC 


,^kk 


“ ^^ko \ “ ^^kko ^ \o ~ \ ‘^ko ^ko ' (4*2.25) 


vjhere 


Tt +r- n, 
12 k 


E-u = 2 

i=l 


j 


li ‘ 




= ^k V ^ 

2 

where = (r^+r^) e^^ . We assume s are independently 

2 

and normally distributed with mean zero and variance e^ . From 
the above assumptions follows a normal distribution with 


= ‘1 - 2Pko ■k’' - <k \o • 

V(Y^o) = (1 - - 2Bi,to ■h '^ko ■ 


(4*2-26) 


In the case of IC Y follows a bivariate normal distribution with 
E(Y) = (I-H)“^ P O ’ 

— -o (i.2.27) 

v(y) = (I-H)“^ (I-H')“^ • 

The extension to a plot containing rows of first crop and 
rows of second crop is straightforward and is on similar lines 
as discussed in section 3*2* 
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4.3 ESTIMATION OF THE PARAMETERS 

(insider the experiment discussed in section 3*2 ^ere 

variation in plant densities has heen introduced by changing 

plant to plant distances within rows* Let d, # d«_ # 

inij^ 2in2 

®1 * 82 be the plant density 

levels in MC for the two species# whereas in IC they are 

# d 


Urn 


2 m ' =’1#2#-.# s* The total number of es^erimental 


plots are KCs^ + S 2 + st) » where K is the number of replications- 
From the equations (4.2.23a/b) individual plot yields of the 
component crops in IC can be expressed by the following 
bivariate sets of simultaneous stochastic equations • 


n - 2Pio rli hlo * ?!! 


"li 

-1 


* ""12 

+ 2 ^ 2 i ^120 '^ 2 m '^ 2 imj ^limj ' 


10 

( 4»3.1a) 


(1 - 2p 


^2 

20 ^2m 


—1 ^22 
^2i ^220 ‘^2m ^ ^2imj 


■ 2 i 

-1 


^21 ‘^2in ’'20 


'21 


^^li ^210 ^limj ^2iinj ' 


C 4* 3 'lb) 


i =1#2#»-# ti m = l#2#..# s# j s=l#2« 


k . 


This can be written as 

-imj ^i °m -imj 


(4.3*2) 
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**21o9 L j-OI" all ths dDssrvations in MC and IC is given by 


-2log L = constant -2K E S log 1 1 - H i + s s log b, „ 

j m 


+ K(s^+st) log V^Q + KCs^+st; log 


20 


+ S S E (G. Y. . - (p)' v“^- (G, Y. . - (p ) 

i m j --o oim im — xmj ^ 

,-l ,,-l Y, . - d. <P, / 




( 4*3«7) 


where 


"im 


5“1 D 1 
i ra 


G._ = P7" D_/- (I-H, _)# V_. = V_ ?.•" D 


dm" ' oim 


D-1 
o i m 


Assuming g_ _ " ^120 ' ^20 ' ^220 ^ ^210^ 

known and equating to zero the partial differentials with 

respect to <Pj^q » ^20 ' ^10 ^*20 ' estimating equations 

of these parameters are the following * 

^ %m ^ ^li '^Im^ ^10 = ^ ^ ^^llim ^limj ■*■ ^12im ^2imj' 
1 m 1 m J 

+ S ' (4.3.8a) 

2 1 I-* 

I P2i ^2m^ ^20 “ ^ ^ ^ ^^22im ^2imj ^lllmYlimj^ 


K( E d^^ + 

m2 2 


+ S E ^20m j ' 

m., J 2 2 


(4-3*8b) 


K(s, +st)V,^ = S E S Pq^^ d^^ ^'°llim "^iimj '*’ ^2im ^2imj ^li^m *^10^ 

^ 1 m 7 


X 2 

<P,«) 


( 4»3*9a) 
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and 


A 

K(s2+st.) V^Q 


-1 .-1 


i m j ^^22im ^2imj + ^21im ^limj“P2l '^2m ^20^ 

(4-3. 9b) 


^ ^^20m2 ^20m2j " '^2m2 *^20^^ ' 


where 

^kim = ^ - - 2(r^^-l) , 

^2im = “2^21 ^120 ^2m^ ’ (4.3*10) 

^21 im “^^li ^210 ' 

Substituting these values in the log - likelihood function 
( 4 . 3 . 7 ) and simplifying / leads to the function 


fCg.) 


constant -2K[ J E log U-H, I +11: log 1+ , ] 

i m k j 

+ KCsj^+st) log V^Q + K(s 2 +st) log V^q * (4.3*11) 


where 


2. = ^^10 ^ ^110 ' ^120 ' ^20 ' ^220 ' ^210^ 

/ 

= , P 2 ' ^3 ' ^4 ' ^5 ' ^6 ^ * 

The ML estimate of g. is the value corresponding to the minimum 
of f(g_) . 
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4*4 FITTING THE MODEL WITH BLOGiC EFFECTS 

By considering the arrangement of the plots in randomized 
blocks * the equations (4«3»la/k)) get modified as • 

Hlim HtaJ = Pll ‘»io + + H2im ^21mj + HlmJ 2 

( 4»4*la) 

^22im ^2imj “ P2i '^2m^‘^20+®2 + ^21im ^limj + ^2imJ * 

(4*4*lb) 

In the case of MC 

^konvj^ ~ ' C4*4»2) 

where a, , ^)^om i effect of the jth block on crop k. 

The log - likelihood function is same as ( 4 * 3 •7) except that 
\o replaced by <P^^. > where <P)coj = ^ ^ko + “^kj^ ' 

eiqjressions in (4»3»8a/b) for estimation of now estimate 

+ 5 / where = 2 a^j/K. 

The es'tima'ting equations for and V 2 Q are 

K(s,+st) V Q = 2 S S ^^llim "^limj + ^12im ^2imj 

i m j 

1 


( 4»4»3a) 



^21 im 


K(32+st) 


= S E S p;^ (b. 


i ;; j -21 '"2m ^“22im ^2iraj 


P2i '^2m *»j)" + 




(4. 4.3b) 


where< 


'*' 1 m ^ ^ ^^llim ^limj'^ ^12im '^2imj^ 




( 4 . 4 , 4 a) 


*^2in "•■ ^ ^ Po-i ^'jTv -i =22 (b 


IT1< 


+ b.,-, Y, ,_,) 


2 i m 21 "2ra' "20 j ' ^ ^ '"221m -21mj + "21im ^limj 


+ 


V 




-20^3 


( 4 . 4 . 4 b) 


The function to be minimized for est.imating g_ is same as 

/\ A 

( 4 . 3 . 11 ) vd.th and V^q as in (4..4.3a/b). 


4.5 MONTE CARLO SIMULATION STUDY 

A simulation study is carried out to examine the performance 
of ML estimators * in the extended model including block parameters. 
The Monte Carlo experiment consists of both IC and MC plots* In 
IC three levels of each factor are selected. They are * 

1:1 ^ 1:2 / 2:1 ; d^jjj = l#2/3 and = 1*2#3. All the 27 
combinations of the levels of these factors are considered . 

In MC there are three density levels for each crop and they are 
d, = 1/2/3 and d^ = 1/2/3 • Randomized block design with 



76 


three replications is used. Hence the experiment consists 
of 99 plots arranged in three complete blocks. 

The values selected for the intraspecific competition 
coefficients are 

^10 = ^20 = ' ^110 = ^220 = “O'SO, = 60.00,^20 = 40.00. 

In IC the interspecific competition coefficients are selected 
as 


^120 ~ ^210 


- 0.10 . 


Values for the error variances and V 
that CV is 15% in MG. We assume c^ = 
^2 ~ ^22 “ ^12 * values considered 

are 

Cl = c^£ = c^^ = ^2 ^ ^^22 = *^12 " ' 


2 Q are chosen such 

Cii = C 2 ^ and 

for these constants 


In each case of =1 and C|,= 2 data for 500 experiments 
are generated by assuming block effects as zeroes* The 
parameters are estimated by using the iterative equations 
given in appendix 4-1 • The initial value are obtained by 
Nelder and Mead simplex algorithm. The convergence criterion 
adopted is same as in section 3*5 . The values of bias 
( 3 . 5 . 1 ) and '^MSE ( 3 * 5 . 3 ) are given in Table 4*1 • For the 
competition coefficients # ' *^jL0 ^20 * relative 

biases do not exceed 2*5% in both the cases* The relative 
to its population values do not exceed 27% when 
c^ = C 2 = l.and 14% when € 3 ^= 0 ., = 2* In general the 
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pararaeters are estimated precisely. However . the size 
of the experiment is considerably large. One can also 
see the precision of the estimates by their asymptotic 
variance-covariances- For example / the asymptotic 
variance — covariance matrix for the competition coefficients 
Var (|.) * when c^ = c^ = 2 , is given below *. 




15 

• 62 9.16 

5. 

72 

-0 .03 

-0.09 

-0.05 





17.42 

6. 

76 

-0.05 

-0 -15 

-0 .08 






6. 

64 

-0 .0 6 

-0 .17 

-0 .09 


Var(i) 

= 10 ' 














19.37 

10.32 

4.30 i 

! 

i 








19.53 

4-25 

I 

1 









1.S7J 

/ { 4 * 5 * 1 ) 

where 

g. = 

‘ho ' 

^110 ' ^1 

20 


20 ' ^ 

220 ' ^ 




^10 

== ^20 

= - 0 . 40 / p 

lie 

> “ 

^220 

-0 .30 

' ^120 "" 

®210 = 


^10 

= 60.00 / <P20 = 

40 . 

00 

' '"10 

= 81.00 

' ’^20 = 

3 6.00 . 


4.6 A PARTICULAR CASE. OF INTEREST 

4 . 6.1 Model in terms of plant densities only 

Investigation of row arrangements and plant density- 
levels simultaneously requires an experiment of large 
size. In many IC experimental pro'^roinmes studies on the 
effect of variation due to plant density levels • at a 
fixed row arrangements^ are not uncommon* A mociel to analyse 
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tne exp0 rimental data from this type of experiment can 
be easily derived from (4*2-22a^n;. Under the assumption 

" ^21 ^2 ^2 ^ ^12 " ^ 


(1 - - 2(ri-l)r-^ 2r;" 0,^ Y, 


4- S, 




.^2 


(1 20 20 ^^2 “ 2 ( r 2 - l ) r 2 ^220 ^2 ^ ^2 ~ '^ 2^20 ^ 21 o 7 l ^2 


-1 


( 4« 6 • lb) 


These equations can be written as 

Cn O , 


(1 - p' dj_h Yj = dj^ »! + p; 


12 ^2 ^2 ^1' 

c. 


( 4-* G * 2s) 


(1 - 02 <^2^^ ^2 = ^2 ^^2 ■*■ ^21 '^l^ ^1 + ^2 ' U.6-2b) 


where 


0^ = 20 


10 


4- 2(r. -I)r7^ 0 


^110 ' ^2 “ ^®20 ^2 ^220 


a' _ 2 t”^ R # <P, = n-, <P- / £, = S S u. ,• ,• ' 

Pmk ~ Pmko k * k xo k 


1 J 


k t= 1/ 2 / i = 1/ 2/ . • l^v ^ 1 = 1' 2i 


. . » n 


k * 


In the matrix form the equations can be written as < 
(I-H) Y = D i + £ ' 


( 4»6.3) 


where 



BO 


(I-H) = 


1 - Pi 




8 ' 

^12 °2 ' 


^ - ^2 ^2 


f = . Y = (Y^, 2^' 


M 2 


The mean and variance are 

E(y) = (I-H) ^ D^ and V(YJ = V (I-H')“^ # 

(4.6.4) 

where -V = diag (V^/V„) / V = p, v . 

X K K XO 

4*6.2 Estimation of parameters 

Let us consider an e^eriment with s density levels 
of both the crops arranged in K randomized blocks* The 
individual plot yields can be represented as / 
c^ C — 

^ *, r 


(1 - d 3 _l) Y^. = d 3 _^ <P 3 _ + d^. Y^i- + e^. .,(4*6. 5a 

c ^ 

(1 - djr Yj. = <Pj + dj_l Y^J_J 


P„ -id + E2^.,C4.6.5b) 


In the matrix notation these equations are / 


(I 

-h> iij = h 2 + 

e . . » 

-IJ 

where 

1 - e( hi 

“Pi 2 ^2i 

(I-H^) 

, , h 

-P'21 hi 

, *"2 ! 
^"®2 ^2i 



' and e j = (e. 


(4.6.6) 
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The equation (4.6.5)* after incluaiiii^ rhe block parameters* 
can be written as 


(I-H,) + qj . 

where <P. = (<P + ^) * a. = (a, .*q; )' 

_j - -j -j Ij 2j 

The value of -21og L is 


( 4 . 6 . 7 ) 


-21og L = constant - 2K E log II-H.I 4 -K s loq(v, .■/„} 

i ^ ^ 


+ SS Id-1 Vj 


(4.6*5/ 


Following the same procedure as earlier * we obcain 

C 

\ = (K + S E E 1(1 - dj^l) Y^.. - S', dj Yj^j! 


1 J 


= (K . E d^i)-! I E 1(1 - S '2 42 ) - S'l d^l Y^ij » 


1^„(K.s)-1eE dj^ , 


k = 1* 2 


where 


®lij 


•2ij 


/S 

hj 


A 

<P . 
2 j 


= Cl - Pi d7) Y,,J - Pi2 d^l Yj^j - d^i . 

, Cl - P'j dip Y^.j - P'21 dl\ Yii^ - d,i J,. 

c ^ 'y 

= ^ - ^I'^li^ -lij ’ ^12 '’ll "Hi' 


(E d^i) E 


{ (1 - P '2 dip Yji, - 6'21 djl Yiij) 

i i 

be minimized for estimation of 3's is 


The function to 
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f^(g_) = log {S S } + log { 


r ^ d_i . , } -2s"^ 2 log D. / 

i j 21 2lj ^ 


where 


i- ^ / -» 2^ 


= 1I-H,1 = (1 - p' d^",)(i -g' d^^) 


"=1 =2 
21 '^ll ^21 • 


(4-6.9) 


4.7 DISCUSSION 

The model suggested by Wright (l-2-3a/b) for plant density 
levels can be obtained by replacing by and by / 
respectively / in the right hand side of the equations 
( 4 . 6 » 2 a/b) / as follows • 

(1 - Pi - Ph 4^' Y = + t ' 


(1 - Pi dh . ah; Vj = Y Y + ^ 


Now taking = '^2 “ ^ obtains 

E(Y^) = d^ - 6' 2 <3^)"^ 

E(Y2) = <32 ~ ^2 " ^21 V ^ ' 

which are same as that of Wright. However, there is a difference 
in the error structure- The assumptions made in extending 
the model for IC to incorporate the variations due to plant 
densities in this chapter are logical extensions of MC 

situation* 

For keeping diseasslons general different oonetants 
=1 ' Yl ' Y2 ' =2 ■ '^22 °21 IntiDduoed in the model. 
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But / in practice/ the constants associated with the density 
may not differ- Hence we can assume c^ = c^^ = c^^ and 

^2 '^22 '^12 ’ situation (4-2.25) c, and c, , 

K KK 

usually lie between 1 and 2* Consequently c, and c.* , are 

taken as 1 and 2 in simulation study. In practice the 
problem of estimation of c^^ and C 2 also arises- It has 
been observed that due to high correlations among the 
estimates of c^ / and competition parameters there is a 
considerable problem in the convergence of the iterative 
process. 

The model developed in section 4*2 involves six 
competition parameters* Use of a good numerical procedure 
becomes very important when the dimension of competition 
parameters is large. It is noticed in some cases that the 
iterating procedure using Newton -Raph son equations converge 
to seme local maximum / which is far off from the values 
given by the simplex procedure. In such cases quasi-Newton 
iterative procedure# which reduces the step length # is 
required to achieve the convergence near the values given 
by the simplex procedure. This procedure was applied to 
obtain the ^^L estimates vhen such a problem arose. 

In general there is a good agreement between the 
estimates of variances obtained oy asymptotic expressions 
(appendix 4-2) and those obtained empirically through simulation 
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Fo r 

example 

^ when 

the parameter 

values are 

^10 

— q — 

- H20 ~ 

-0 . 40 / 

^110 = 

P 2 20 “ -0 . 30 4 

®120 ~ ^210 

o 

1 — 1 
& 

= 60.00 

II 

o 

40.00 , 

o 

1! 

00 

• 

O 

o 

'/ = 4:0 « 00 


~ ^^2 “ 

■che empirical estimates and estimates using asymptotic 
expressions are the following I 

Estimated Variances of the t^aroineter x 10"" 



o 

1 — 1 

oa 

hio 

^120 

® 20 

^220 

3 

210 

Empirical 

0*22 

0 .17 

0.0 2 

C . 27 

0*15 

0 

.01 

Asymptotic 

0.16 

0.17 

0 -07 

0.19 

0 . 20 

0 

.02 


From the simulation results * and also from the 
asymptotic variance - covariances # it is evident that 
the parameters can be estimated precisely for the type of 
experiment considered in section 4*5 at 15^ uV . oow’ever# 
the sizes of the experiment and blocks are considerably large 
It is/ therefore / important to search for certain optimum 
designs such that the parameters can be estimated precisely 
through an experiment of reasonable size* 
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NUMERICAL METHOD FOR FINDING THE t-lAXIMUM LIKELIHOOD SSTII-LLTES 


For simplicity / taking = S 2 the ML estimate of 
g_ correspond to the minimum of the function 


-1 .-1 _2 


f^B) = log ( B E E 


1 .2 




i m j 




+ I E E p;l =1, + E^ E 


-2( S S log D^^ + S 2 log ) / g 


i m 


k rn., 


kom 


( A 4 • 1 ♦ 1 


k 


where 


£ (g_) = (f(g_) - constant) K ^ g ^ > 
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APPENDIX 4.2 

variances - covariances 0? THE ESTIMATORS 


The expiressions for the asymptotic variance - covariance 
matrix for the estimators of the parameters are derived 
in the following ^ 

The log - likelihood of the observations from IC and 
MC is 
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i m ^ ^ 
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where 
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value to each crop. One can think of an Index such as 

Monetary Equivalent Ratio (MER) , „hlch takes into account 

the monetary value of IC and also of the corresponding MC- 

LER and MER can be defined in tontis of expected yields and 

denoted by LER and MER / as 
e 0 

.laER^ =s {ECy^)/E(Y^^)) , (5.1.1a) 

K K 

MER^ ={l\ E(Y^) } /{ S R^^ P E(Y_; } , (5.1.1b) 

where R, is the cost per unit Y, / and o, = r^ ( r-, 4 .r..) . 

A. X K ‘"12 

LER^ and MER^ will be approximately equal to E(LEr) and 
E(MER) • MERg not only depends on the expected yields but 
also on the price ratio R 2 ./R 2 * When R^/R., = S( Y.^^ ) ,/Et Y^^) ^ 
MER^ is identical to tER^. MER^ represents the monetary 
advantage in IC compared to the situation in which the 
component crops are grown in MC in the same proportion as 
in IC . Instead of MER^ one can also use the variate total 
monetary value (TMV) / given by 

TMV = S { 

}c 

Calculation of bER^^ and MER^j_ for the IC with r^^ * 
row arrangement based on the model (2.2*8) is straight 
forward . When the effect of plant density is also incorporated 
the problem of selecting optimum MC yields in calcularing 
these indices also arises- For the kth MC we have from (4-2. 2-5 

= =k e 

^^^ko^ = "^ko ‘^ko^^ - Sk \o ^ 


(5.1-2) 
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where ^kko^' “^k the family of the 

curve* 


lim E(Y,_^^) = 0 fo 

■ko" ° 


y.o- =k 


lim eCy^^) = - 0y^ »hen = 1 . 

‘ko* ” 


= o 


when > 1 


Hence the curve given by (5*1*2) is asymptotic when c^ - 1 

and ECY, ) attains maximum yield given by 
KIO 

Whereas in case of c^ > 1 . attains maximum at a finite 

-1 /c^ 


d^Q ' say # given by dj^ = ^ 

n (j* will be negative and is not of interest in our 
^k 

case. Hence when > 1 the optimum yield of monocrop k is 

obtained at the density d^^ • When c^ = 1 » tne optimum 

yield is given by (- • These can be used in calculating 

the indices hER^ / MER^ etc* 

5. 2 METHODS AND EXAMPLES OF OPTIMIZATION 

5.2*1 Optimum ,rgw arrangement 

5. 2.1.1 Methods of findiW row arranqement 

4 - r- : r , of the component crops 
For the ith row arrangement u 2i 

in IC / we have from (2.2.8) » 

E(Y^) = ^ Pi (i • 

s -1 . 

For the kth MC eCY^^) = ^''^Pkk 
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^ OT" r,^ • "*0 or / the IC svstem tends to MC- From 
'hen ■^]_3_^- 2 i 

■he above expected yields / the values for are / 

'eli 


= E(Yj^j^) / ElYj^g) 

= (1-2^11) I l- 2 (r 2 -lJ ■(• 2 r;i Pj^ Pj ] 

( 5 . 1 . 3 a) 


I-e2i = ■-' 


. ( 1 - 2 S 22 ) d ‘^[ ! l-2(ri^-l)r];i hi ’^21 + Pli ®21 *"1 ^ 2 ^ ^ 

( 5 • 1 . 3b) 

= II-H.l = (l-2p3_3_+2r"i 03_l)(l-2022+2^2i®22^"2r^^r2^B3_2623 

Let us consider certain cases of practical interest 

, -1 „ 

and ^22= ^12 ^mk 


where 


i) When = P 21 

we have / 


^ 6„v 'k m = 1-2. 


Cl-2B^l) dT^ [^1-2^22^ + 2P22Cr^itr2iHr^.T2^ 


)(r. Pli 


'eli 


(5.1. 4a) 
. -1- 


I, „ = Cl-2h,) Pih[(l-2hl) P 2hi(PlVP21>‘Pll-P2i'‘ ] Pai ' 

(5.1.4b) 


where 


Dj 


^ (l- 23 il) (1-2022^ 2rii3iiU-2t522 


(1-20.,2^ 2r2i ^22 


-1 6..C1-23 ii)- 


T J- Hi ♦ i. f O XT 

From this it follows that > eli e2i 

e.erY i - api thus there is no biological abvanta.e of 


over MC for eny row 


arrangement. 
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ii) When 8 


11 ~^22~^12“®21' then 


L 


eki 


Pki 


every i , i-e- the expected yield of any component 
crop in IC / irrespective of row arrangements » is 
same as the expected yield to be obtained through MC 

iii)when = 822 ' then 


.-1 


ler^i = oyr (1-25,2 ((1-20,1 + ^diht’Pil 


-1 

‘li ^^11' ■^2i 


+ 11 - 23j^^ 4 - 2r-,^ 0-nJ - 1 - 


"li^^2i' 


(012 ^2 + ^21 ^ 


wnere 


p, , (1 - 25,, 4. 2r"i p„)(l - 2i3„ + 2r;^ B,,) 


-1 


2( r2_ j* r2j_) ^21' 

From the above relation it follows that LER value remains 

e 

same when and r2£ sre interchanged for any i* 

The coefficients ^ ^12 * ^22 ^21 

degree of intra and interspecific competitions- Apart from 
the cases mentioned above LER^ depends uoon the row arrangements 
in general- The row arrangements which are of interest in 
practice are very limited as ana - 2 x take only integer 


values • Moreover / 


ind ^ increase^the component 


crops in IC tend to behave like MC- One can easily aOtk 

out LER . » MER^. / TMV. or any other a' E(Y^) for any i- 
Comparison of these values will help in arriving at. the 
optimam strategy about the row arrangements- ?or any fixed 
proportion -there are many possible row arranger,* 


lents* Eor 



example . when . = p^^ = 1/2 , ^he k 7 possible row 
arrangements are given by r^^^ = = 1 , 2, . . . , n . instead 

of considering variation in . and r 2 ^ a better insight 

could be had by studying the effects due to the variation 
in geometry and proportion p^ . v/e shall denote the effect 

of geometry by g / where g is the highest common integer 

among such that (r^^ : r^^) = 

example when = 1/2 / for ; r 2 j_ = 1 • 1 / g^ = 1 / 


for r^Li : ^21 = ^ 


• 2 / = 2 / when Pij^ = 2/3 / for 


^li • ^2i “ ^ 


3 . g. = 1 . for r^. ; r^^ = 4 : 6 , g. = 2- 


As is the propoxition of the expected hth crop yield 

in. IC for the ith row arrangement to its expected .'-IC yield / 
^^eki ^ ^ki^ indicates the degree of advantage for crop k. 


'eli 


> Pli leads to the condition 


^^12 “ ^ 11 ^ - f ^® 11^12 “ ^ 11^22 ^21 ^^ 11^22 " ^ 12 ^ 21 ^^' 

A similar relationship holds for i>^ 2 i * 


5 • 2 • 1 . 2 Some examples 

Example 1. Let the values of the parameters be = ^22 “ 

3 i 2 = 0 21 = !^-i = ' 5 '^ ' ^2 ^ ' 

The values of LER^ and MER^ worked out from (S-l-la.b) for 
different row arrangements are given in Table 5.1a- In this 
case interspecific competition is sufficiently weak in 
comparison to intraspecific competition for both the crops- 
LER^ is maximum for 1 '• 1 arrangement as is evident from 

It remains same even if we interchange r 3 _ and r 2 


Table 5.1a • 



Table 

5*lb Values of tER 

and MKK 

V JLH pet 




^1 

^21 

at the 

— —0 * 40 / 1 

parameter 

= 60.00 

values 

/ h., =. 

3 2^ • 40 

50.00 . 

' ^12 = 

-O . lO t 8 

^2 

1 

2 

3 

4 

6 

8 

10 






, 

, , - 


1 

0*945 

1*068 

1.063 

1.055 

1.042 

1*033 

1.028 


(0*936) 

(1*063) 

( 1 . 0 60 ) 

(1*052) 

( 1 . 0 40 ) 

(1.032) 

(1.025) 

2 

1*078 

1*081* 

1.068 

1.057 

1.044 

1.035 

1 .0 30 


(1*075) 

(1.079) * 

(1.066) 

(1.056) 

(1.043) 

(1.034) 

(1.028) 

3 

1*075 

1*069 

1.059 

1*051 

1*040 

1.033 

1.028 


(1*073) 

(1.068) 

(1*058) 

(1*050) 

(1 .040) 

(1.032) 

(1.027) 

4 

1*065 

1*060 

1.052 

1*046 

1*037 

1.031 

1.026 


(1*064) 

(1*059) 

(1;051) 

(1*045) 

(1*036) 

(1.030) 

(1.026) 

6 

1*050 

1*046 

1*0 41 

1.037 

1*031 

1.027 

1 .023 


(1*050) 

(1 .046) 

(1.041) 

(1*037) 

(1.031) 

(1.026) 

(1.023) 

8 

1*040 

1*037 

1.034 

1.031 

1.027 

1.024 

1.021 


(1*040) 

( 1 *037) 

(1 .034^ 

(1.031) 

(1*027) 

(1.023) 

(1.0 21; 

10 

1*034 

1.031 

1*029 

1.0 27 

1.024 

1.0 21 

1*019 


(1*033) 

(1*031) 

(1.029) 

(1*027) 

(1.023) 

(1.021) 

(1.019) 


Note* * indicates the maximum values of and . 
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because = ^22 * However » partial LER values L, 

e 1 e 

or ^2e difi.erent when rj,^ and r^ are interchanged- The 

parameters / 3 j _2 ' ^22 ® 21 '^Haracterise the degree 

of competition- For this case p'^ = - 0.0 83 and 821 = - 0 - 120 . 
Though 3 2 ^ ~ ^22 1^12 ~ ^ 21 * depend upon rj^ 

and r 2 since p3_2 ^21 * Table S-la MER^ has been 

calculated for = 2 and = 3* In this case MER^ is also 
maximum for 1 • 1 arrangement ♦ Reason for such a conclusion 
is not difficult to see from the values of the coefficients 
011 » 022 ' ^12 ^21 suggest that the competition 

pressure per unit area is less when the neighbouring rows 
belong to the other component crop than of its own type- 


Example 2- Let the values of tne parcimeters be 0^^. = * 

^ -0-60 / 0^2 = -O-IO' 021 = 60 / = 50- 

For this set both LER^ and KER^ (Table 5 - lb) are minimum 
fQ^ 1*1 row arrangement and are less tnan unity i-e- MC 
is advantageous over 1 • 1 row arrangement- In tnis case 2-2 
arrangement is better than the other arrangements- Here 
0^2 = 0*083 and 021 = 0 - 480 . The value of 02 i is considerably 
large in comparison to its value in example 1- Tnis is 
reflected in LER^ and MER^ values for various row arrangements. 

From Tables S-la^b it is clear that as r, and r 2 increase . 
LER^ and HER^ tend to unity i-e. they tend to the values of 
MC L doth the exan.ples. Enures 5.1 and 5.2 give the values 
of LER^ and HER^ for different values of geometry (gi of the 


different 
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component crop rows for the two sets of parameters- 

indicates the degree of intimacy between two component 

crop rows and it varies between zero and one. Intimacy 

is maximum when g ^ is one and minimuiti vhen it is zero* 

From figuresf.l and 5.2 it is evident that LER and HER 

e e 

tend to unity as g increases- These figures give insight 
into the dependence of bER^ and on the intimacy 

between the two component crops- bER^ and MER^ vary more 
steeply at the lower values of g but rhe changes at the 
higher level of g are small* For the first set of parameters 
the values of bER^ are symmetrical around p^ =0.5 / 
consequently figures are drawn only for p^^ < 0.50. 


5.2*2 Optimum row arrangement and plant density 

Let us now consider the problem of finding simultaneously 
the optimum combination of row arrangement and plant density 
levels for both the component crops which maximize Z = E(Y) . 
Z may be LER » TMV or any other meaningful linear combination 
of the component crop yields- Since d^ and d 2 take continuous 
values and Z is a complicated function in d^ and d 2 ^ we are 
required to use some nonlinear optimization methods for 
finding optimum dj_ and d 2 ,‘ When = ^ 2 ^ and 

= 2 = =22 = =12 ' 


E(Y) = (I-H) D 9 # 


where 




l-^.d 


^2 - 


2 2 


^21'^ 


^ 12^2 ^ 

^ # Y = (Y^#Y2) * 

1-81%^ , t' " 



LERe 



p^=0-75 Bl^=O-40,J3i2--O-lO,J32f-O-6O,J32f-O-40 

~D jJl^ 60-00, /J2= 50-00 

1-06 

-Pl=0-10 ^ 

1-02 - 

0-98- 

0-94 -p^-0-50 

0-90^ ' ^ ' i ' 6 ' 8 ' lb 

Geometry (g) 

Fig.51 Values of LERe 'or different geometry (g) ond 
proportion p^ 
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f3^^ = -0-40,13i2=-0-^0.n22='0‘60 

J32i = -0-40,jJi = 60-00) jJ2=50-00 


t 

10 


0-98- 

0-94 

0-90 


p^ = 0-50 


0 


4 


1 

6 


8 


L ^ 

Geomet ry (g) 

F,g.S2 values ot MERe tor dlHcreni 

- proportion p, at the prices R,=Z and Rj 
crop 1 and crop 2 respectively 
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are examined graphically, when 8,^-6 o ^ 

f12 - P 21 = ' Cl = = 1 

and Qj^ and d^ tend Infinity , z approaches asymptotioally 
to the maximum yield -(a^ S'l + <e^ p-^) . it tends 

to zero v.'hen and C2 > 1 * Nature of the function is such 
that when 82^2 P 21 close to zero the optimum and 

d2 occur at larger values of and 02 • This is also 
reflected in Tables 5 . 2 a,b and c- When c^ = 1 optimum 
occured at relatively larger value of d^ compared to 
the situation when c^^ > !• In practice the parameters 
have to be estimated from the experimental data* In situations 
•when the optimum lies outside one raiiye of the experiment it 
is advisable to repeat the experiment by including levels 
around the predicted optimum and d^- 

The optimum row arrangement for Zj_ / Z2 and LER^ can be 
found by comparing their values over different r^^ and r2 
obtained by maximizing over and d2* Tor the linear 
combination and LER^ the optimum row arrangement is 1*1 ^ 
V7hile for Z2 1*3 arrangement is optimum* For example » the 
nature of the curves of » Z2 and LER^ for variation in 
dj, and d2 at their optimum row arrangements / are shown in 

Figures 5.3 (c^ = “"s = * '=2 = 

from Figure 5.x that Zy 23 and LER^ do not attain their 

p. .3 m anh n - 10 » whereas in Figure 5*3 

maxima even when d^ = 10 ana 02 - 

. r..r..o 4 ...ca relatively smaller values of 

maximum value is attained at rciarivej.y 


d^ and d2‘ 






Table 5* 2a Maximum of Z 


/ Z 2 and LER^ at their optimum densities d * and d^ , 
for various row arrangements- 0 io=p 2 O=“O* 4 O '3 -R 

-^.01, = 60.00 , <p^ , 40.00, o,“°c' Jf . 


*^120 “ ^210 


10 


= ~o • 30 / 


Row arr 


1 

1 

1 

2 

2 

2 

3 

3 

3 

4 


1 

2 

3 

1 

2 

1 

2 

3 

4 


ts 

^1 



^2 



LER 

e 


* 

^2 

2l 


^ * 
^2 

z * 

^2 


a * 
^2 

LER*' 

e 

1.08 

1.13 

38. 5S^* 

0.97 

1.12 

93-63 

1.13 

1-14 

1 . 30 ' 

1.06 

0.96 

36. 28 

0 . 0 

0.96 

99. 27 

1.12 

0.96 

1.18 

1.06 

0.92 

35 . 53 

0 . 9'^ 

0.92 

103 . 57 ** 

1 .07 

0*94 

1.13 

0.94 

1.11 

33.71 

0.91 

1.12 

70,58 

0.97 

1.06 

1.18 

0.94 

0.95 

33 . 20 

0 . 69 

0.96 

80 . 57 

0.89 

0.99 

1*12 

0.94 

0.91 

33.22 

0.91 

0.91 

87.75 

0.95 

0.95 

1.09 

0.90 

0 . 93 

31.56 

0 . B9 

1.1'^ 

59 . 22 

0.90 

1.03 

1.13 

0.91 

0.96 

31.56 

0. 8B 

C . 9 

69. 48 

0.92 

0.92 

1.09 

0.92 

0.93 

31.84 

0 . SB 

0.91 

77.31 

0.93 

0.88 

1-08 

0.89 

0.89 

31. 23 

0.85 

0.89 

75.82 

0.83 

0.92 

1 • 0 6 




Note. ** indicates overall optimum VrAue 


Ztl + 2EU,,) 

'■^1 " "“^ l ^ *^ l ' '^ 2 ^ ' '^2 


^2 - S^Y^) + 7E(Y2) 

ler* = l.ER^(a*.dp 
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. ' ' ' ''"fV., i ■■’' ■ ' , ’ ’■ 


In this .chapter some aspects of optimizing certain 
functions of component crop yields such as > MER^ 


n of the component crop 


^4 if. j In practice it may be 

I j ' 4 - ‘ ' '''* 

ife^ejirfeain constraints which component 


4 'a^g,’ested by Mead and Sterm (1979) * 
11984) may be more appropriate* These 


can be worked out in terms of the expected 
case of LERg and MSRg. 

Use of nonlinear optimization methods in finding the 


optimum plant density levels require some care* One may 


have to verify that the maximum thus obtained is a local 
maximum or a global maximum* The functions considered 
here are generally concave functions of d^ and d 2 and 
there exists a unique maximum * This has also been verified, 
by plotting the values over different d^ and d 2 * In practice 
the cost of inputs are also to be considered hence# optimization 
subject to these economic constraints may be ^ mo re appropi 


The numerical results along vdth the graphical repres.e^^'^bhs 
vjill help in understanding the response of component 

. 'V,' ' I ^ I ''‘ii 

yields to various row arrangements- and density levels * This 
will subsequently help in arriving at their optimum values* 
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In practice the parameters have to be estimated through 
experimental data« The pb^cision of the estimates depends 
on the amount of information available in the data* Large 
sample variance ~co variance matrix and simulation results 









ANALYSIS OF EXPERIMENTAL DATA 


6.1 INTRODUCTION 


fitted adei|uately by smaller number of parameters* For example# 

R sen %he, component crops behave independently in IC# the 
usefulness of interspecific competition coefficients has to 
be tested. Likelihood rati© criterion^and test based on 
asymptotic standard errors^can be used for testing whether 
the comoetition parameters are different from zero or not* 


Akaike (1974) information criterion (AIC) 


number of parameters to be fitted in the model. This < 
adjusts the log-likelihood for the number of pararaete 
be estimated. AIC is derived under the assumption tha 
true distribution can be described by the given model 


Table 5 

. 2c Maximum of 

Z. / Zt and LER^ at 
12 c 

their optimum densities d* and d* for 

va irious 



arrangements 

• PlO = ^20 = 

' ^110 = 

^220 = 

P 120 

= ^210 = 

- 0 . 01 , 



*^10 

= 60.00 

/ <p^ = 40.00 e 

1 and 

£= 2 • 


Row 

arr 

_ ■ 

Z-, 


z 




- 

angements 

1 


^2 



ler 

e 


^1 

^2 

^1 

d* 

2 



* 

"2 

* 

d* 

2 

ler* 

e 

1 

1 

7.86 

1.38^*^ 

50. io 3*36 

1. 21 

99.35 

6.15 

1.39 

1. 28 

1 

2 

7. 34 

1.08 

43*84 3. 23 

1.01 

103.00 

6.18 

1.00 

1.16 

1 

3 

6.18 

0.96 

41.11 ?.1B 

0.94 

106. 34** 

6.14 

0.95 

1.12 

2 

1 

9.91 

1 . 38 

45.26 4.59 

1.22 

77.85 

3.95 

1 . 40 

1.13 

2 

2 

9. 20 

1.03 

41.73 4 . 40 

1.00 

85.94 

9.02 

1.03 

1.08 

2 

3 

9.13 

0.99 

40 .01 4.33 

0.95 

93.02 

6.51 

0.83 

1.05 


1 

12.67 

1. 47 

43-77 5.55 

1.24 

67.72 

10. 61 

1.51 

1.08 

3 

2 

11.11 

1.08 

41.23 5.36 

1.01 

76. 20 

10. 29 

1.00 

1.04 

3 

3 

1C. 7 2 

1.00 

39.87 5.31 

0.95 

82. 38 

10. 30 

1.00 

1.03 

4 

4 

12* 25 

0.96 

39-09 6.05 

0.92 

81. 60 

10.33 

1.00 

1.02 

Uo te 

• ** indicates overall optimum values 
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5 . 3 DISCUSSION 

In this chapter some aspects of optimizing certain 

functions of component crop yields such as LER # MER 

e e 

or any other linear combination of the component crop 
yields have been discussed. In practice it may be 
desirable to impose certain constraints which component 
crop yields must satisfy • In sucn situation modified 
Indices of LER suggested by Mead and Sterm (1979) / 

Reddy and Chetty (1984) may be more apompriate. These 
modified indices can be worked out in terms of the expected 
yields as in the case of LER and MER . 

Use of nonlinear optimization methods in finding the 
optimum plant density levels require some care. One may 
have to verify that the maximum thus obtained is a local 
maximum or a global maximum. The functions considered 
here are generally concave functions of dj^ and d^ and 
there exists a unique maximum . This has also been verified 
by plotting the values over different d^^ and d^* In practice 
the cost of inputs are also to be considered hence# optimization 
subject to these economic constraint may be more appropriate- 
The numerical results along with the graphical representations 
•will help in understanding the response of component crop 
yields to various row arrangements <a.nd density levels • This 
will subsequently help in arriving at their optimum values* 
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In practice the parameters have to be estimated through 
experimental data* The precision of the estimates depends 
on the amount of information available in the data* Large 
sample variance-covariance matrix and simulation results 
discussed in the third and fourth chapters give some idea 
on these aspects* Extrapolation should be avoided as the 
model discussed here may not hold good for very high values 
of plant density and row arrangements- 



CHAPTER VI 


analysis of experimental data 

6.1 INTRO EUCTION 

Eventhough the model developed in the second and fourth 
chapters are based on the principles of plant competition the 
goodness of fit of these models to the experimental data 
collected under field conditions need to be examined. As 
the models involve a large number of parameters# it is of 
considerable importance to examine whether the model can be 
fitted adequately by smaller number of parameters* For example' 
v/hen the component crops behave independently in ICt the 
usefulness of interspecific competition coefficients has to 
be tested* Likelihood ratio criterion^and test based on 
asymptotic standard errors^can be useu for testing whether 
the competition parameters are different from zero or not* 

Alternatively# Akaike (1974) information criterion (AIC/ # 
defined by 

AIC = -2 log L + 2(number of parameters) # 
appears to be a useful technique in judging the adequate 
number of parameters to be fitted in the model. This criterion 
adjusts the log-likelihood for the number of parameters to 
be estimated. AIC is derived under the assumption that the 
true distribution can be described by the given model when 
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its parameters are suitably adjusted. The model with the 
minimum AIC gives the best fit. Examination of residuals 
is another way of looking at the goodness of fit of the 
model. Unlike the Likelihood ratio or Akaike measure this 
does not depend upon the form of the alternative model. 

In section 6.2 some aspects of fitting the model developed 
in the fourth chapter which include the factors such as 
row arrangement and plant densities shall be discussed. The 
data are from an IC e> 5 )eriment on mustard and chickpea 

(Kushwaha/ 1983 i • The experiment involves a very limited 
number of row arrangements and plant density levels. These 
data have been used for testing the goodness of fit of the 
model developed in previous chapters# in the absence of data 
from better designed experiments* 

6 . 2 ANALYSIS 

5.2.1 Experimental detail s_ 

The e:cperiment has sixteen treatments with four replications 
each. Randomized block design was adopted- Among tne sixteen 
treatments four treatments are of MC on mustard# four treatments 
are of MC on chickpea and the remaining eighc are on IC of 
the above two* In each MC there are four plant density levels* 
In IC there are tvXD row arrangements and four density levels 
for each row arrangement. In MC plots interrow distances 
are 60 and 30 cms for mustard and chickpea* respectively. 
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The tvjo row arrangements in IC are : r^^^ =1 -4/2*2- 

where r^^^^ and r 2 ^ are the number of rows of mustard and 
Ohickpea in the ith row arrangement- In IC plots the 
distance between a row of chickpea and mustard is 45 cm. 

The area occupied by one row of mustard is equal to the 
area occupied by two rows of chickpea- The row arrangements 
in MC and IC are illustrated in Figure o-l. The plant 
density levels are obtained by adjusting the space between 
the plants within a row- The four olant density levels 
for mustard and chickpea are : 

2 

p, = P_ z= 7.5/ 15*0/ 22-5 ana 3a. O plants/m 
Im 2m Jr- / 

They have been coded as d^j^ and d^^^^ given by 


= a. 


= 1/ 2/ 3/ 4 


where 


d-,. 


Im ^ "2m 

= P- /7.5 / k = 1/ 2. 

xm knr 

The plant density levels for each crop in MC are same as 
in IC . Plotwise yield data are given in appendix 6-1. 


6.2*2 Pitting of ^the model 

The data have been analysed by fitting the models ( 4 . 3 . 1 a/b) 
and (4.3*4a/b) for IC and HC / respectively • The models are 
as the following * 
for IC 


^^“^im^ -imj “ ^i °im -^j -imj 


i =1/2 / m/j =l/2/.-/4 
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for MC 


(1 ~ “ ^^kko ^konij,j “ *^koj ^kom^j 


k = 1,2 ; m^ . J = 1,2,3, 4 ; 


where H 


CP . , E . 


, / Uj ,F.,,-_. 

im im i — oj — imj — imj 


'%ai 


are 


as mentioned earlier* Here we assume Cj^ = and 

*^2 ~ ^2 2 ~ ^12 * present experimental situation 

'^Im = ^2m = ^Im^ = '^ 2 m 2 m , m^ , = 1,2, 3, 4 ,* 


-1 


Pli " ^^li^^*li ■*■ ^21^ ' P2i " ^2i^^^li ^ ^2i^ 


-1 


The log -likelihood for all the data is given by 


log L =-{48 log(2t ) + 2( E log d^^ + S log d^^) 


m 


m 


2m' 


+ 2 S S logCp^^ P2i ^im ^2m^ ^ 
i m 


+ 4 £ S log 1 1-H .^1 + 4 r 1 
1 m km 


- 24 logCv^Q.V^o) 


I £ £ £ [ ! (I-H,„) 2oj !' 

1 m j 


. {(I-H. ) Y. . - P. bo ]] 

xm — imj i ra — oj 


I E E £ [ { - d^ 

k m j 


ko j 


r ] 


where 


b. 


km 


1 + f*VVo' 


>d:o' Ik 


The parameters are estimated by using the iterative 
procedure discussed in appendix 4*1* The initial values 
are obtained by simplex procedure. The convergence 
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criterion used is same as in section 3*5- Iterative procedure was 
repeated for various values of c^ and since there was a problem 
in the convergence for these data in obtaining the I'iL estimates of 

and c^ • The value of -2 log L is near its minimum when c^ = c^ = 1.5- 

The estimates of the parameters with their asymptotic standard 
errors are : 


Parameter 6^^ g 


Estimate -0.428 -0.17B -0.069 C.07S -0.526 -0.036 3.142 6.370 0.223 0*407 

Asymptotic 

standard 

error 0.100 0*101 0.023 0.046 0.065 0*278 0*353 0*40'^ C.054 0-094 
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A A 

The correlations are very high among -r 
among 0^0 ' ^220 ^ 


and also 


’210 


whi cn ociii 


oe seen from the following 


correlation matrix Corr*(p) (given in the upper triangle) 

1.00 


Corr. (g) 


-0.67 

“0 . 2 

-O.IB 

0.16 

0.19 

1.00 

0,91 

0 .21 

-O.IB 

-0.22 


1.00 

0.21 

-0.19 

-0.22 



1.00 

-0.88 

-0.94 




1.00 

0 .84 





1.00 


The adequacy of tne competition parameters is also examined 
by dropping one parameter • two parameters etc.# at a time by 
likelihood ratio criterion- The values of -2 log L and AIC under 
different constraints on the parameters are ; 





Cons 

traints 


1 

— - i 

is 

o n s z r o. — , 
int |S 

1 ' ' ‘ 

1 

10=^110=^120= j 
20=®210=^’210'^| 

‘ l'>0^ ! 

1 

' ' t ' ’ ' “ ? 

‘ i g 

^20 ty ^ ^ 

j , "'21C* 

20^ 

^210=^' 

-2 log L 

1 ; 

j 198.909 

! 

530-361 ! 

A ^ i 

Cui. • o o ' 

! 

200 . 87lil99.40 6 1204-15 2 

201.971 

AIC 

j 230.909 

550-361 

! t 

1231 .5 '3 6 ! 

230. B71 1 2 29 .40 5 {23 2.15 2 

’■) og ,071 


The value of -2 log L decreased considerably by fitzing all 
the competition parameters* AIC is minim’om for the model fitted 
with the constraint ^ 71 q =0* 
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Hence this model is considered for furtner analysis though 
the increase in AIC is small when 3^20 ^210 

taken as zeroes. The estimates of the parameters when 62 x 0 ' 
with their asymptotic standard erixjrs are 


Parameter 


^10 


'110 


120 


20 


^220 10 


20 


'^lO 


20 


Estimate -O. 4 I 6 -0-192 -0-077 O-O 4 O -43-491 3-149 6-825 0-219 0-432 


Asymptotic 

standard 

error 02.00 0-102 0-024 O-Olo O.O 36 0- 352 C- 30 G 


-C96 


The observed and expected yields, averaged over the repli^^ 

with the residuals for both the crops in AC and IC,are 1 

The goodness of fit of this model is examined by 

residuals against: me expected values (Picure 6 * 2 ;* Prom 

• « • • » 1*1 '*“h. '*"he 

figure it is evident tnat tne variability is increasing '-* ' 

expected yields# which is eicpected from the model. Otne^^^® 

residuals do not show an’/ pattern and thus the model app^ 

satisfactory for these data- In some of mu=^tard the oa uo ao 

not follow smooth pattern as expected / this is perhaps drought 

experienced by this crop (Kushwaha/ 1983; - Consequently ' some o- ^he 

residuals are large. 



0'40- 


175 


0 - 20 - 


Mustard 


0 - 0 ' 


0 - 20 - 


• • t 

• « 

•• !♦ • 

•0*60 


• voo 


1-50 

Expected 
, yields (t /he 


0-60- 


C hick pea 


)'40- 


0 - 20 - 


0 


1-00 


0 - 20 - 


• • 


2-00 






! \ • 

3-00 


• •] 

4-00 

• Expected 

,* , yieldslt/ho) 


0-40' 


0-60- 


Fig.6-2 Residuals against the expected yields 
tor the fitted model 



Table 6«1 Observed 


’lant 

iensibies 

1 

2 

3 

4 

1 

2 

4 

1 

2 

4 


Observed 

1 .614 
1-233 
1-233 
1-274 

0.639 
0.514 
0.526 
O . 464 

0.95S 
0.784 
O .80 2 
0.771 
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diiQ 0xp>0c*t.Go msGn vicX'i<5 3n,<'5 ir0sxo.u.3ls I t/hs.} • 


Mustard 


Chickpea 



Expected Residual 

s Observed 

Expected 

Residual s 

1.422 

MC 

0.192 

3 . 464 

3-587 

-O .123 

1.419 

-O .186 

3*723 

3.717 

0 .00 6 

1 .292 

-0 .059 

3-717 

3 -574 

0 -143 

1.197 

0.077 

3 . 466 

3-321 

0 -145 


IC 





1*4 




0.516 

0 .123 

2.876 

2-746 

0 . 1‘30 

0.5 23 

-0 .009 

2-971 

3.18 4 

-0 .213 

0 . 47 6 

0.050 

3-032 

3-093 

—0 -0 61 

0.432 

0 .032 

2-917 

2.909 

0.008 


2^2 




0.976 

-0 .OIP 

1.753 

1. 612 

C - 1 41 

0 .961 


2-341 

2.103 

u • 3 3 

0 .860 

-0 .0 58 

2-216 

2-425 

-0 . 20S 

0.772 

-0.001 

2.023 

2-121 

-O -093 
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6.3 DISCUSSION 

Significance tests for the estimates based on their 
asymptotic variances and likelihood ratios give similar- 
conclusions* However' Akaike s criterion appeals to be a 
useful tool for selecting the parsimonious models although 
it does not indicate whether the better of the two models 
is significantly better. Problems in convergence for 
obtaining the ML estimates iteratively have been experienced 
in certain cases* This is because or nigh correlations 
among the estimates due to limited number of row arrangements 
and density levels used in the experiment. Hence one has to 
look for certain designs such that the correlations among 
the estimates is small. Detaile'"^ analysis of residuals has 

not been attempted because of tne limited data. Looking at 

vxV '■■/'/.i 6.:' 

the plot of the residuals there ^Tiijure 6.2.* does not appear 
to be any pabbern and bne fit appears sati s facto ry. 

r * ^ 1 

However/ to establish the goodness of rit or tne models 
developed here more ereperiments need to be aricilysed. 


The estimates of the competi tion coetficients/ broadlv 
indicate the nature and degree or competition among the 
crops* The interspecific competition is considerably smaller 
than intraspecif ic competition among tne prams of both the 
crops* The estimates of intraspsciric competition coefficients 
are 


h ” hio^ 


-0.603' 1 2 = ^^20 ^ h20' 


= -0.451 
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for mustard and chickpea/ respectively. Competition among 
mustard plants is more than in chickpea since 3^^ = I .3532 • 
The estimates of interspecific competition coefficients are 


6 ' = <P 

P120 20 


-1 

(f 


10 ^120 -^* 166 / 3 


210 


e>- ^ » 

^ <P P - O i 
10 20 ^210 ~ * 


for mustard and chicJcpea / respectively . There exists a 
significant interspecific competition for mustard whereas 
it is negligible for chickpea. 


Optimum combination of row arrangement and plant densities 
can be obtained by adopting the methods discussed in the fifth 
chapter* In the case of MC optimum density levels for obtaining 
the maximum yields are 

d*Q = 1*35 4 and d^ = 1*700 / 

for mustard and chickpea/ respectively. The maximum ercpccted 
yields corresponding to these densities are 1*47 and 3*87 f/ha* 
Here *320 is positive which is rather unexpected. This is 
perhaps due -to limited amount of information available on 
the row arrangements. However / this is nearly nonsic'niticant - 
In extrapolation of the yield corresponding to row arrangements/ 
particularly when r^ =1/ P 20 a dominant role. When 

it is positive then the expected yield for chickpea/ when 
r 2 = 1 behaves in a very unexpected fashion/ and consequently 
it has not been considered for optimization purposes. 
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The optimum densities in IC which maximize (5.1*la) 
for various row arrangements # except when r^ = 1 / are 
the following ’ 


Table 6.2 Optimum densities d* (mustard) and d*( chiciqsea) 


. No • 

^1 

^2 

^1 


* 

LER 

e 

1 

1 

2 

3-601 

2 .043 

1.326 

2 

1 


2. 283 

1.595 

1.235 

3 

1 

4 

2.276 

1-660 

1.190 

4 

2 

2 

2 .00 3 

1.609 

1.197 

5 

2 

3 

1.987 

1.666 

1.153 

6 

2 

4 

1.980 

1.718 

1.128 

7 

3 

2 

1.867 

1.636 

1.143 

8 

3 

3 

1.854 

1.687 

1.115 

9 

3 

4 

1 .848 

1.736 

1.099 

10 

4 

2 

1.788 

1.648 

1.112 

11 

4 

o 

1.776 

1 .697 

1.092 

12 

4 

4 

1.772 

1.744 

1 -081 


These results give some indication for the future 
direction of experimente>.tion* as wt, increase r^ and r2 
decreases- Among the above combinations 1*2 row arrangement 
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the range of the observed densities/ and they depend upon 
the row arrangement. They are relatively larger when the 
crops are more intimate/ consequen-cly the values of LER^ 
are also larger* This indicates that more investigations 
including ill arrangement are required- 



Plot yields (t/ha) of each crop over the row arrangements .. r 2 and plant densities 
P, f (plants/m^) (Kushwaha/ 1983 ) » 
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Notation 


^kk 

^km 

^ko 

^kko 

^kmo 


autoregressive coefficient represents the 
effect of crop k row on the same crop 
k m = 1,2 autoregressive coefficient represents 
the competition effect of crop m on crop k 
autoregressive coefficient represents the competition 
effect between two neighbouring plants of crop k. 
represents the competition between the two neighbouring 
rows of crop k at plant density of crop k dj^ = 1 
k 7 ^ m = 1,2 represents the competition of crop m on 
crop k when the plant density of crop rn^d^ = 1. 


Corrections 


1 . Page 9 and 10 

2. Page 13, line 
just above 
(1.2.13) 

3. Page 21 line 3 
line 7 

4. Page 23 line 6 


5. Page 24, (2.2.9) : 

6. Page 25,(2.2.11) : 

line 9 i 

line 15 : 


7. Page 26 

line 14 


the section numbers should be 

1.2. 2. 3 and 1.2. 2. 4 instead of 2. 2. 2. 3 and 

2. 2. 2. 4. 

o 2 

erf not a 


Vshen * n ot * whe re 

^ni,rj^+r2+l ^m,r+r2+l 

The correct specification of Yj^ and are 


n 


ri+r 

Y. = E 
^ i=l 

yii» 

Y 2 = 2 y2i $ 

2 i=rj^+l 

n 


ri+r2 

e = 2 

uii , 

’ ^2 = 2 U2^ 

1 1=1 

i=rj^+l 


Hh’ 1 not H*A1 


“ki “li 


u,.^ not s 


"ki 


N. 


ki 




I'i. 


u. 


® 1 = 2 "li * 
^ i=l 


N^+N2 

^2 = ^ , ^2i 

^ i=N +1 ^ 

2 r , “2i 

i=i’* +i 


where = brj^, k = 1,2. 


In the matrix B, (1,2) element is p,, not 


^111 

^ , n and u not y , n and u . 
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